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BACKGROUND 


Residual  stresses  are  those  contained  in  a  body  which  has  no 
external  traction  or  other  sources  of  stress,  such  as  thermal  gradients 
or  body  forces.  When  the  body  is  externally  loaded,  these  stresses  are 
called  internal  stresses,  and,  accordingly,  residual  stresses  may  be 
considered  as  a  special  case  for  vanishing  external  loads.  Residual 

stresses  result  from  non-uniform  plastic  deformation  which  includes  cold 
working,  forming,  forging,  heat  treatment,  etc.  Their  presence  in 
manufactured  components  plays  an  important  role  in  determining  the 

behavior  of  the  component  when  it  is  subjected  to  service  loads  and 
environment.  ~Tt  has  also  been  shown  that^1  /he  residual  stress 

distribution  directly  affects  the  growth  rate  and  frequency  of  formation 

of  stress  induced  cracks  in  steels.  - ___ 

Only  in  the  case  of  surface  stresses  in  components  made  of 

crystalline  materials  can  nondestructive  evaluation  of  stresses  be 
performed  by  the  X-ray  diffraction  method.  Although  considerably 
improved  in  the  last  ten  years,  this  method  still  suffers  from  serious 
problems  which  severely  restrict  its  aoppi ications.  Ultrasonic  methods 
appear  to  hold  the  best  promise  in  measurements  of  bulk  stresses  in  both 
crystalline  and  non-crystalline  materials.  Calculations  have  shown  that 
ultrasonic  velocity  changes  are  linear  functions  of  applied  stress  and 
combinations  of  second-  and  third-order  elastic  constants.  In  the 


application  of  these  calculations  to  determine  unknown  stresses,  both 
the  velocity  in  the  absence  of  stress  as  well  as  third-order  elastic 
constants  have  to  be  known  independently.  In  addition,  the  measured 
velocity  strongly  depends  on  microstructural  features  which  makes  it 
necessary  to  develop  a  calibration  between  velocity  and  stress  in  order 
to  be  used  in  the  determination  of  unknown  stresses.  Development  of 
preferred  orientations  (texture)  during  deformation  or  fatigue,  severely 
modify  the  third-order  elastic  constants.  These  problems  can  be  solved 
when  the  differences  between  velocities  of  shear  waves  polarized 
perpenducular  and  parallel  to  stress  direction  are  used.  Due  to  these 
differences,  a  shift  in  phase  will  occur,  and  the  out-of-phase 
components  will  interfere  and  cause  a  change  in  intensity.  This  method, 
however,  does  not  have  at  present  enough  sensitivity,  and  requires  an 
accurate  determination  of  the  shear  velocity  in  the  absence  of  stress. 

Basically,  the  temperature  dependences  of  the  elastic  constants  are 
due  to  the  anharmonic  nature  of  the  crystal  lattice,  and  can  be  related 
to  the  pressure  dependence  of  these  constants.  If  we  consider  the 
isothermal  bulk  modulus  By  to  be  a  function  of  pressure  and  temperature, 
it  follows  that 

*BT,  _  1  r3BT,  SBT.  .  (1) 

^It  “  jjFj:  gT~l v  Trip1 


where  3  is  the  volume  coefficient  of  thermal  expansion.  Swenson  has 
shown  empirically  that  for  many  materials 


and  it  can  be  shown  by  differentiating  (B§  -  By)  that 
3B-  3Bt 

IpIt  “  TpIt 


(3) 


to  an  accuracy  of  a  few  percent.  Thus  it  follows  that  if  Swenson's  rule 
is  correct, 


3BS  -1  3BT 
pH  T  =  pT  *  pH  P 


(*) 


to  an  accuracy  of  a  few  percent.  The  right-hand  side  can  be  calculated 
from  the  measurements  of  the  temperature  dependence  of  the  second-order 
elastic  constants.  The  left-hand  side  is  calculated  from  the  third- 
order  elastic  constants  C.jk,  and  eqn.  (3)  is  then  expressed  as, 


i! 

3T 


I  1C111  +  3C112  +  3C113  *  2C123’ 


(5) 


More  rigorous  relationships  can  be  derived  for  the  temperature 
dependences  of  the  longitudinal  and  the  shear  moduli  which  are  related 
to  the  temperature  dependence  of  the  ultrasonic  velocity  v  as, 


3ft  nM  _  nV 

TT  rTT 


(6) 


From  this  general  argument,  it  can  be  seen  that  the  temperature 
dependence  of  the  ultrasonic  velocity  (longitudinal  or  shear)  is  a 
measure  of  the  anharmonic  effects  generated  when  the  solid  is  subjected 
to  a  stress.  Changes  in  the  anharmonic  properties  due  to  the  presence 
of  residual  stresses  can  therefore  be  detected  by  the  changes  in  the 
slope  of  the  relationship  between  the  ultrasonic  velocity  and  the 
temperature.  This  slope  can  be  determined  with  a  high  accuracy  as  its 
value  depends  only  on  measurements  of  the  relative  changes  in  the 
velocity  and  not  on  the  absolute  values.  Experiments  undertaken  on 
aluminum,  copper  and  A533B  steel  elastically  deformed  showed  that  the 
ultrasonic  velocity,  in  the  vicinity  of  room  temperature,  changed 
linearly  with  temperature,  and  the  slope  of  the  linear  relationship 
changed  considerbly  as  the  amount  of  applied  stress  was  varied.  In 
aluminum,  copper  and  steel,  the  relative  changes  of  the  temperature 
dependence  of  longitudinal  velocity  increased  respectively  by  24%  at  a 
stress  of  96  MPak,  by  6$  at  a  stress  of  1 80  MPa,  and  by  15%  at  a  stress 
of  240  MPa.  The  results  also  showed  that  the  temperature  dependence  of 
ultrasonic  velocity  (dV/dT)^  at  applied  stress  o  may  be  represented  by 
the  relationship, 

(dV/dT)  -  (dV/dT)o 

- R573 TT0 - 


where  (dV/dT)o  is  the  temperature  dependence  at  zero  applied  stress  and 
K  is  a  constant.  The  accuracy  of  the  temperature  dependence  in 


measuring  applied  stress  is  estimated  to  be  +18  MPa  in  aluminum,  +25  MPa 
in  copper  and  MPa  in  steel. 

All  the  above  studies  have  been  performed  when  the  stress  is 

applied  in  a  direction  perpendicular  to  that  of  the  ultrasonic  wave 

propagation.  In  practice,  the  magnitude  as  well  as  the  direction  of 

applied  or  residual  stress  are  not  known.  In  order  to  apply  the 

temperature  dependence  of  ultrasonic  velocity  method  to  stress 
measurements,  relationship  similar  to  that  expressed  by  eqn.  (7)  should 
be  available  for  stress  applied  parallel  as  well  as  perpendicular  to 
wave  propagation.  It  is  also  important  for  the  application  of  the 
method  that  the  theoretical  background  for  these  relationships  is 
established  so  that  values  of  the  constant  K  in  eqn.  (7)  can  be 
calculated  using  the  solid  anharmonic  constants. 


SUMMARY  OF  RESULTS  AND  DISCUSSION 
~r - - - - - - 

The  research  on  this  contract  was  initiated  in  Juiy  1982  in  order 
to  develop  the  method  of  the  temperature  dependence  of  ultrasonic 
velocity  for  the  nondestructive  evaluation  of  bulk  residual  stresses. 
Results  previously  obtained  on  aluminum,  copper  and  steel  have  shown 
that  this  method  offers  a  promising  possibility  for  measurements  of  bulk 
stresses  in  these  materials.  The  results  also  indicate  that  the  method 
is  less  sensitive  to  variations  caused  by  alloy  composition  and  perhaps 
texture.  Three  objectives  have  been  set  for  the  research  performed 
during  the  two  year  period  of  this  program.  The  first  objective  is  to 
establish  experimental  relationships  between  the  temperature  dependence 
of  ultrasonic  velocity  and  stress  applied  in  a  direction  parallel  to  and 
perpendicular  to  that  of  wave  propagation.  These  relationships 
represent  the  variations  of  the  temperature  dependence  when  the  stresses 
are  applied  in  tension  as  well  as  in  compression. 

Appendix  I  includes  two  papers  which  describe  the  results  obtained 
on  the  effects  of  applied  stress  on  the  temperature  dependence  of 
ultrasonic  longitudinal  velocity  in  aluminum  6061-T6.  The  first  paper 
published  in  the  Proceedings  of  the  1982  IEEE  Ultrasonic  Symposium, 
deals  with  the  study  performed  to  establish  the  relationship  between 
ultrasonic  velocity  and  compressive  stress  applied  perpendicular  to  wave 
propagation.  The  results  show  that  the  temperature  dependence  decreases 
linearly  by  as  much  as  20%  when  a  stress  of  80  MPa  is  applied  in  this 
direction.  The  results  also  indicate  that  the  relative  changes  in  the 
temperature  dependence  in  aluminum  6061-T6  are  equal  to  those  previously 


obtained  on  other  three  aluminum  alloys,  and  empohasize  the 
insensitivity  of  the  temperature  dependence  method  to  alloy  composition. 


The  second  paper  in  Appendix  I  which  appears  in  the  Proceedings  of 
the  14th  Symposium  on  NDE  (1983),  discusses  the  effects  of  stress 
applied  in  a  direction  parallel  to  and  perpendicular  to  ultrasonic  wave 
propagation  on  the  temperature  dependence  of  longitudinal  velocity  in 
aluminum  6061-T6.  The  results  obtained  in  this  study  show  that  the 
temperature  dependence  of  ultrasonic  velocity  increases  linearly  with 
either  tensile  or  compressive  stress  when  they  are  applied  in  a 

direction  parallel  to  that  of  ultrasonic  propagation.  In  the  case  of 
stress  applied  perpendicular  to  propagation,  the  results  Indicate  that 
the  temperature  dependence  decreases  linearly  with  both  stresses.  The 
calibration  curves  constructed  to  relate  the  relative  change  in  the 
temperature  dependence  to  applied  stress,  give  an  estimate  of  ±10  MPa 
for  the  sensitivity  of  these  curves  in  determining  unknown  applied 
stresses  in  the  aluminum  used. 

The  second  objective  of  the  program  was  to  employ  previously 
determined  relationships  between  temperature  dependence  of  ultrasonic 

velocities  and  applied  stress,  to  measure  residual  stresses  of 
predetermined  values.  Appendix  II  includes  a  preprint  of  a  paper 
published  in  the  proceedings  of  the  ASM  Symposium  on  Applications  and 
Development  of  Nondestructive  Evaluation  for  Use  in  Materials  Processing 
held  in  Philadelphia,  October  1 983 »  in  which  measurements  of  the 
temperature  dependence  was  used  to  determine  the  stress  distributions  in 
a  disc  made  of  A533B  steel.  The  stresses  in  the  disc  is  introduced  by 

shrink  fitting  a  rod  of  the  same  steel  in  a  hole  drilled  in  the  disc 


with  a  Slightly  smaller  size.  The  temperature  dependences  were 
determined  using  longitudinal  and  shear  waves  propagating  along  the 
thickness  but  polarized  either  parallel  to  or  perpendicular  to  a  radial 
direction.  The  results  indicate  that  the  temperature  dependence  of  the 
longitudinal  velocity  remains  unchanged  as  a  function  of  the  radial 
distance  on  the  disc,  which  agrees  with  predictions  of  elasticity  theory 
for  the  behavior  of  the  axial  stress  component.  On  the  other  hand,  the 
temperature  dependences  of  the  shear  velocities  show  considerable 
variations  with  the  radial  distance,  reflecting  the  changes  in  the 
radial  and  the  tangential  components  of  stress  as  a  function  of  the 
radial  distance. 

The  third  objective  of  the  program  was  to  establish  the  theoretical 
background  for  the  relationship  between  the  temperature  dependence  of 
ultrasonic  velocity  and  stress.  Experimental  results  obtained  on  three 
materials  (aluminum,  copper  and  steel)  show  that  this  relationship  is 
linear  and  to  a  large  extent  does  not  depend  on  the  alloying  composition 
in  each  material.  The  experimental  results  also  show  that  the  slope  of 
this  relationship  is  large  in  aluminum  and  small  in  copper.  Appendix 
III  includes  a  preprint  of  a  paper  to  be  published  in  the  proceedings  of 
198**  Review  of  Progress  in  Quantitative  NDE,  which  describes  the 
derivations  of  the  temperature  dependence  of  ultrasonic  velocity  with 
respect  to  stress,  and  the  comparison  between  theory  and  experiment. 
The  derivations  are  based  on  the  expressions  obtained  by  Hiki,  Thomas 
and  Granato  for  the  temperature  dependences  of  isothermal  elastic 
constants  of  single  crystals  using  a  quasiharmonic  anisotropic  continuum 
model.  The  temperature  dependence  of  ultrasonic  velocity  in  a 


polycrystal  is  dervied  in  terms  of  second-,  third-  and  fourth-  order 
elastic  constants  using  the  Voight  average  of  Young's  modulus  and 
rigidity  modulus.  This  expression  along  with  the  experimental  values  of 
the  temperature  dependence  of  longitudinal  velocity  obtained  in  aluminum 
and  copper,  yield  values  of  fourth-order  elastic  constants  which  are  in 
good  agreement  with  those  computed  by  other  investigators.  Furthermore, 
the  derivative  of  this  expression  with  respect  to  stress  is  shown  to  be 
constant  for  constant  values  of  the  stress  dependence  of  the  third-  and 
the  fourth-order  elastic  constants.  This  result  agrees  with  the 
experimental  findings  relating  the  temperature  dependence  of  ultrasonic 
velocity  to  stress. 
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Abstract 

The  temperature  dependences  of  ultrasonic  ve¬ 
locities  are  due  to  the  anharmonic  nature  of  the 
crvstal  lattice,  and  therefore  can  be  used  for 
jtress  measurements.  Experiments  performed  on 
t>06l-T6  aluminum  show  that,  in  the  vicinity  of  room 
temperature,  the  ultrasonic  longitudinal  velocity 
decreases  linearly  with  temperature,  and  the  slope 
0f  the  linear  relationship  varies  considerably  when 
the  specimen  is  subjected  to  stress.  For  compres¬ 
sive  stress  applied  perpendicular  to  wave  propaga¬ 
tion,  the  temperature  dependence  of  the  velocity  is 
found  to  decrease  linearly  by  as  much  as  20%  at  a 
stress  of  80  MPa.  The  results  also  indicate  that 
the  relative  changes  in  the  temperature  dependence 
of  velocity  as  a  function  of  stress  are  equal  to 
those  previously  obtained  on  other  aluminum  alloys. 
This  shows  the  insensitivity  of  the  temperature  de¬ 
pendence  method  to  texture  and  alloy  composition. 
The  method  thus  offers  a  promising  possibility  for 
the  nondestructive  measurement  of  residual  stress. 


1.  Introduction 

There  is  a  general  agreement  that  ultrasonic 
methods  appear  to  hold  the  best  promise  in  the  non¬ 
destructive  measurements  of  bulk  stresses  in  both 
crystalline  and  non-crystalline  materials. * Cal¬ 
culations  have  shown  that  changes  in  ultrasonic 
velocities  are  linear  functions  of  applied  stress 
and  unknown  stresses  can  be  determined  when  both 
the  velocity  in  the  absence  of  stress  as  well  as 
third-order  elastic  constants  are  known  indepen¬ 
dently.  The  measured  velocity,  however,  strongly 
depends  on  microstructural  features  which  makes  it 
necessary  to  develop  a  calibration  between  velocity 
and  stress  for  each  material  in  order  to  be  used  in 
the  determination  of  unknown  stresses.  In  addition, 
development  of  preferred  orientations  (texture) 
during  deformation  or  fatigue,  severely  modify  the 
third-order  elastic  constants  to  be  used  for  the 
calibration.  These  problems  can  be  solved  when  the 
differences  between  velocities  of  shear  waves 
polarized  perpendicular  and  parallel  to  stress 
directions  are  used.  Due  to  these  differences,  a 
shift  in  phase  will  occur,  and  the  out -of -phase 
components  will  interfere  and  cause  a  change  in 
intensity.  This  method,  however,  does  not  have  at 


present  enough  sensitivity,  and  requires  an  accu¬ 
rate  determination  of  the  shear  velocity  in  the 
absence  of  stress. 

Basically,  the  temperature  dependences  of  the 
elastic  constants  of  a  solid  are  due  to  the  an- 
harmonic  nature  of  the  crystal  lattice. 3, 4  a 
measure  of  the  temperature  dependence  of  the  ultra¬ 
sonic  velocity  can,  therefore,  be  used  to  evaluate 
bulk  stresses.  Experiments  undertaken  on  aluminum 
and  copper^, 6  elastically  deformed  in  compression 
showed  that  the  ultrasonic  velocity,  in  the  vicin¬ 
ity  of  room  temperature,  changed  linearly  with 
temperature,  and  the  slope  of  the  linear  relation¬ 
ship  changed  considerably  as  the  amount  of  applied 
stress  was  varied.  In  aluminum,  the  relative 
changes  of  the  temperature  dependence  of  longitu¬ 
dinal  velocity  decreased  linearly  by  as  much  as 
23%  at  a  stress  of  approximately  96  MPa.  The 
results  obtained  on  different  types  of  aluminum 
alloys  also  indicate  that  the  relative  changes  in 
the  temperature  dependence  due  to  applied  stress 
are  insensitive  to  composition  and  texture,  and 
the  data  obtained  on  these  alloys  can  be  repre¬ 
sented  by  a  single  relationship.7 

All  the  above  studies  have  been  performed 
when  the  stress  is  applied  in  a  direction  per¬ 
pendicular  to  that  of  the  ultrasonic  wave  propa¬ 
gation.  In  practice,  the  magnitude  as  well  as 
the  direction  of  applied  or  residual  stress  are 
not  known.  In  order  to  apply  the  temperature  de¬ 
pendence  of  ultrasonic  velocity  and  stress  should 
be  available  for  stress  applied  parallel  as  well 
as  perpendicular  to  wave  propagation. 

In  this  paper  the  effects  of  applied  compres¬ 
sive  stress  on  the  absolute  as  well  as  the  tem¬ 
perature  dependence  of  the  ultrasonic  longitudinal 
velocity  have  been  investigated.  The  experiments 
were  performed  with  the  stress  applied  in  a  direc¬ 
tion  parallel  to  and  perpendicular  to  the  ultra¬ 
sonic  wave  propagation. 

2.  Experimental  Procedure 

The  specimens  used  in  the  present  study  were 
made  from  one  inch  diameter  bar  stock  of  commer¬ 
cial  6061 -T6  aluminum  in  the  form  shown  in  Fig.  1. 
All  specimens  were  made  identically  except  the 
overall  length  L  was  varied.  The  specimens  were 
machined  with  a  2.54  cm  diameter  cap  on  each  end 
which  allowed  the  same  specimen  to  be  used  for 
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Fig.  1  Specimen  used  for  ultrasonic  velocity 
measurements . 


stress  applied  in  compression  as  well  as  in  tension. 
The  two  caps  at  the  ends  were  made  flat  and  parallel 
to  within  ±0.002  cm  in  order  to  avoid  diffraction 
effects  in  the  ultrasonic  beam  during  propagation. 
These  caps  were  also  connected  to  the  center  por¬ 
tion  of  the  specimen  by  a  0.06  cm  radius  in  order 
to  minimize  stress  concentrations.  After  experi¬ 
ments  of  the  stress  applied  parallel  to  the  axis 
of  specimen  were  completed,  two  parallel  surfaces 
of  thickness,  t,  were  milled  in  the  center  of  the 
specimen  to  allow  measurements  for  stress  applied 
perpendicular  to  wave  propagation  direction.  In 
this  case  the  ultrasonic  waves  were  propagated  in 
the  center  portion  of  the  specimen  where  the  stress 
is  expected  to  be  uniform. 

The  application  of  external  stress  was  carried 
out  with  a  model  112S  floor  type  Instron  machine 


Fig.  2  Loading  system  used  for  the  application  of 
compressive  stress  in  a  direction  parallel 
to  the  ultrasonic  propagation. 


Fig.  3.  Loading  system  used  for  the  application  of 
compressive  stress  in  a  direction  perpen¬ 
dicular  to  the  ultrasonic  propagation. 


The  application  of  external  stress  was  carried 
out  with  a  model  1125  floor  type  Instron  machine 
of  20,000  kg  maximum  load  capacity.  Two  different 
types  of  loading  arrangements  are  used  in  the  pre¬ 
sent  work.  Figures  2  and  3  show  the  systems  used 
for  applying  the  compressive  stresses.  To  help  en¬ 
sure  the  uniformity  of  stress  in  the  specimens, 
special  effort  was  made  in  designing  these  systems 
to  minimize  effects  of  misalignment  between  the 
axis  of  the  specimen  and  the  loading  axis.  A 
linear  ball  bearing  served  as  the  first  alignment 
between  the  upper  and  lower  loading  axes  and  a 
hemispheric  steel  ball  served  as  the  second  align* 
ment  between  the  specimen  and  the  loading  axis. 

The  shafts  used  in  the  compression  tests  were  mad# 
of  surface  hardened  steel  rod  (Rc  =  40)  to  resist 
possible  wearing  from  the  bearing  during  loading. 

The  ultrasonic  velocity  was  measured  using  the 
pulse-echo  overlap  method  which  has  been  fully  “•* 
scribed  elsewhere. 8  Figure  4  displays  the  expert* 
mental  system  used  in  this  investigation.  A  Pul** 
of  approximately  1-y  sec  duration  of  variable  p**1 
repetition  rate  is  generated  by  the  ultrasonic 
generator  and  impressed  on  a  transducer  of  a 
mental  frequency  of  10  MHz  which  is  acoustically 
bound  to  the  specimen.  The  reflected  rf  echoes  ^ 
received  by  the  same  transducer,  amplified,  and 
played  on  the  screen  of  an  oscilloscope.  Two  ot 
the  displayed  echoes  are  then  chosen  and  exactly 
overlapped  by  critically  adjusting  the  frequency 
of  the  cw  oscillator.  This  frequency  f,  accura 
determined  by  the  electronic  counter,  is  employ 
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Fig.  4:  Fulse-echo  Overlap  System. 


to  compute  the  ultrasonic  velocity  using  the  rela¬ 
tion  V  =  2if,  where  i  is  the  length  of  the  specimen. 
X-cut  transducers  are  used  for  the  generation  of  the 
longitudinal  waves,  which  were  used  in  all  measure¬ 
ments  . 

The  temperature  control  system  is  designed  to 
enclose  the  specimen  and  its  gripping  assemblies  in 
order  to  ensure  stabilized  temperature  for  the  whole 
specimen  during  the  time  required  foT  velocity 
measurements.  These  measurements  were  performed 
between  300  K  and  370  X  where  the  coupling  condition 
between  the  transducer  and  the  specimen  was  found 
to  be  satisfactory.  At  all  temperatures,  the  actual 
temperature  of  the  specimen  was  measured  by  a  cop- 
per-constantan  thermocouple  attached  directly  to 
the  specimen.  The  thermocouple  along  with  a  poten¬ 
tiometer  provided  a  measurement  of  the  temperature 
with  an  accuracy  of  *0.01  K. 


3.  Results 

Temperature  Dependence  of  Ultrasonic  Velocity: 

The  effects  of  applied  stress  on  the  tempera¬ 
ture  dependence  of  ultrasonic  longitudinal  velocity 
were  investigated  on  three  specimens  1,  2  and  3  of 
6061 -T6  aluminum.  Specimens  1  and  2  were  used  when 
the  stress  was  applied  parallel  to  wave  propagation, 
while  specimens  2  and  3  were  used  when  the  stress 
was  applied  perpendicular  to  the  direction  in  which 
ultrasonic  waves  were  propagated.  Typical  examples 
of  the  results  obtained  on  specimen  1  showing  the 
variations  of  ultrasonic  velocity  with  temperature 
at  the  compressive  stresses  0  and  -19.8  MPa  applied 
parallel  to  the  propagation  direction  are  shown  in 
Fig.  5.  From  this  figure  it  can  be  seen  that,  in 
the  temperature  range  from  310-370  K,  the  ultra¬ 
sonic  velocity  decreases  linearly  with  temperature 
and  the  slope  of  the  linear  relationship  varies  when 
the  stress  is  applied.  A  computer  program  was  used 
to  process  the  velocity-temperature  data  to  deter¬ 
mine  the  temperature  dependence  of  ultrasonic  ve¬ 
locity  (dV/dT). 


Table  1  lists  the  results  of  (dV/dT)  obtained 
on  the  three  specimens  investigated  when  the  stress 
was  applied  parallel  to  (specimens  1  and  2)  as  well 
as  perpendicular  to  (specimens  2  and  3)  the  direc¬ 
tion  of  wave  propagation.  Also  included  in  this 
table  are  the  offset  uncertainties  which  measure 
the  degree  of  linearity  of  the  relationship  between 
velocity  and  temperature.  From  these  results  one 
can  see  that  the  temperature  dependence  of  ultra¬ 
sonic  velocity  in  aluminum  increases  or  decreases 
according  to  whether  the  stress  is  applied  parallel 
to  or  perpendicular  to  the  direction  in  which  the 
ultrasonic  waves  are  propagated.  The  results  also 
show  that  the  changes  in  the  temperature  dependence 
are  larger  when  the  stress  is  perpendicular  to  wave 
propagation  than  when  the  stress  is  applied  parallel 
to  propagation  direction. 

Because  the  values  of  (dV/dT)  at  zero  applied 
stress  were  found  to  vary  among  the  specimens  in¬ 
vestigated,  the  relative  change  in  the  temperature 
dependence.  A,  due  to  the  application  of  stress 
was  calculated  and  its  values  are  listed  in  column 
6  of  Table  (1) .  The  variations  in  the  temperature 
dependence  of  the  three  specimens  tested  at  zero 
stress  are  believed  to  be  due  to  differences  in 
the  residual  stress  in  these  specimens.  The  values 
of  A  were  calculated  using  the  relationship. 


(dV/dT) o  -  (dV/dT) 
(dV/dT) 


where  (dV/dT)  is  the  temperature  dependence  at  zero 
applied  stress. 

The  relative  changes  in  the  temperature  depen¬ 
dence  of  ultrasonic  velocity  as  a  function  of  ap¬ 
plied  stress,  obtained  on  the  three  specimens  in¬ 
vestigated  are  plotted  in  Fig.  6.  The  plots  show 
that  the  data  points  obtained  on  specimens  1  and 
2  when  the  stress  was  applied  parallel  to  wave  pro¬ 
pagation,  lie  on  a  straight  line  which  passes 
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Table  (1)  -  Variations  of  the  temperature  dependence  of  longitudinal  ultrasonic  velocity 
with  applied  compressive  stress  in  aluminum  6061-T6. 


Specimen 

Number 

Specimen  Length 

cm 

Applied  Stress 
(MPa) 

dV/dT 

(m/s-K) 

Offset 

Uncertainty 

(m/s) 

A\ 

STRESS  //  PROPAGATION 

1 

3.897 

0 

-0.969 

.166 

0.0 

-19.8 

-0.982 

.234 

1.5 

-39.7 

-0.996 

.252 

3.0 

-59.6 

-1.016 

.233 

5.0 

-79.5 

-1.032 

.342 

6.7 

2 

6.985 

0 

-1.009 

.226 

0.0 

-9.9 

-0.975 

.272 

1.0 

-19.8 

-0.984 

.230 

2.0 

-39.7 

-0.998 

.289 

3.4 

-S9.6 

-1.022 

.174 

5.8 

STRESS  1  PROPAGATION 

2 

t  =  0.624 

0.0 

-1.007 

.358 

0.0 

L  =  6.985 

-16.6 

-0.966 

.326 

-3.9 

-33.1 

-0.933 

.358 

-7.1 

-66.2 

-0.863 

.197 

-14.1 

3 

t  =  0.624 

0.0 

-0.890 

.332 

0.0 

L  =  7.638 

-22.1 

-0.845 

.234 

-5.3 

-44.1 

-0.779 

.313 

-12.6 

-66.1 

-0.736 

.211 

-17.5 

through  the  origin  with  a  slope  equal  to 
-8.63  x  10"4  per  MPa  and  a  correlation  coefficient 
equal  to  -0.992.  The  figure  also  shows  that  the 
results  obtained  when  the  stress  was  perpendicular 
to  wave  propagation  (specimens  2  and  3),  also  lie 
on  a  straight  line  which  passes  through  the  origin, 
but  with  a  slope  equal  to  24.4  x  10“  per  MPa  and 
a  correlation  coefficient  of  0.982. 


perpendicular  to  wave  propagation.  The  slope  of 
the  straight  line  representing  the  velocity  and 
stress  when  the  ultrasonic  waves  are  propagating 
parallel  to  stress  is  -0.210  m/s -MPa  with  a  cor¬ 
relation  coefficient  equal  to  -0.998.  The  slope, 
however,  is  equal  to  0.0832  m/s-MPa  with  a  cor¬ 
relation  coefficient  equal  to  0.993  when  the 
stress  is  applied  perpendicular  to  wave  propa¬ 
gation. 


Stress  Dependence  of  Ultrasonic  Velocity: 

The  effects  of  applied  stress  on  the  velocity 
of  longitudinal  ultrasonic  waves  propagating 
parallel  to  and  perpendicular  to  stress  has  been 
measured  on  specimens  1  and  3  respectively.  The 
measurements  were  performed  at  room  temperature 
with  applied  stress  ranging  from  0  to  -59. 6  MPa 
for  specimen  1  and  from  0  to  -44.1  MPa  for  speci¬ 
men  3.  The  results  of  these  measurements  are 
given  in  Table  2,  which  also  lists  the  values  of 
the  relative  change  in  the  velocity  AV/V.  Plots 
of  these  values  as  a  function  of  applied  stress 
for  the  two  specimens  are  shown  in  Fig.  7.  From 
this  figure,  one  can  see  that  AV/V  increases  or 
decreases  linearly  with  stress  according  to 
whether  the  stress  is  applied  parallel  to  or 
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Table  (2)  -  Variations  of  longitudinal  ultrasonic  velocity  with  stress  in  specimens  1  and  3  of 
6061 -T6  aluminum  at  room  temperature. 


Stress  //  Prop  (Specimen  1) _ Stress  J.  Prop  (Specimen  3) 


Stress 

(MPa) 

Velocity 

(m/s) 

AV/V 

Stress 

(MPa) 

Velocity 

(m/s) 

(AV/V) 10^ 

0.0 

6118.6 

0.00 

0.0 

6191.2 

0.00 

-9.9 

6119.7 

1.83 

-5.5 

6191.1 

0.18 

-19.9 

6121.8 

5.28 

-11.0 

6190.6 

-1.03 

-29.8 

6123.7 

8.40 

-16.5 

6190.2 

-1.70 

-39.7 

6125.7 

11.6 

-22.1 

6189.5 

-2.78 

-49.7 

6128.1 

15.5 

-27.6 

6188.9 

-3.76 

-S9.6 

6130.7 

19.7 

-33.1 

6188.2 

-4.86 

-38.8 

6187.5 

-6.06 

-44.1 

6186.9 

-7.01 

Stress  Dependence  dV/do 

« -0.210  m/s -MPa 

Stress  Dependence  dV/do 

=  0.0832  m/s -MPa 

Correlation 

Coefficient 

-  -0.998 

Correlation  Coefficient 

«  0.993 

Fig.  S  Effect  of  applied  compressive  stress  on 
the  temperature  dependence  of  ultrasonic 
longitudinal  velocity  in  the  aluminum 
alloy  6061-T6.  The  stress  is  applied 
in  a  direction  parallel  to  the  ultra¬ 
sonic  propagation. 


Fig.  6  Effect  of  applied  compressive  stress  on 
the  relative  change  in  the  temperature 
dependence  of  longitudinal  ultrasonic 
velocity  in  6061 -T6  aluminum.  Stress 
was  applied  parallel  to  and  perpendicular 
to  wave  propagation. 
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Fig.  7  Relative  change  in  ultrasonic  velocity 
as  a  function  of  applied  compressive 
stress  in  6061 -T6  aluminum.  Stress 
was  applied  parallel  to  and  perpen¬ 
dicular  to  ultrasonic  wave  propagation. 


4.  DISCUSSION 

1.  Temperature  Dependence  of  Ultrasonic  Velocity 

The  effects  of  compressive  stresses  on  the 
temperature  dependence  of  longitudinal  ultrasonic 
velocity  have  been  studied  in  the  aluminum  alloys 
2024-0  and  6063-T4  by  Salama  and  Ling®.  These 
experiments  were  performed  with  the  stress  applied 
in  a  direction  perpendicular  to  the  ultrasonic 
propagation.  The  results  obtained  by  these 
authors  show  that  the  relative  change  in  the  tem¬ 
perature  dependence  of  ultrasonic  velocity  is  a 
linear  function  Of  the  applied  elastic  compressive 
stress  and  can  be  represented  by  the  relationship 

(dV/dT)  -  (dV/dT) 

a - (3^5 -  -  Ko  (2) 

where  (dV/dT)  is  the  temperature  dependence  at 
zero  applied  stress,  (dV/dT) 0  is  the  temperature 
dependence  at  an  applied  stress  e,  and  K  is  a  con¬ 
stant  equal  to  2.3  x  10"3  per  MPa. 

In  the  present  investigation,  the  experiments 
were  performed  to  study  the  effects  of  compressive 


Fig.  8  Relative  change  in  the  temperature 
dependence  of  ultrasonic  velocity  as 
a  function  of  applied  stress  in  three 
aluminum  alloys. 

stress  on  the  temperature  dependence  of  longitudinal 
ultrasonic  velocity  on  three  specimens  (1,2  and  J) 
of  the  aluminum  alloy  6061-T6.  The  relative  changes 
in  the  temperature  dependence  of  longitudinal  ultra* 
sonic  velocity  as  a  function  of  applied  stress  ob¬ 
tained  in  this  work  as  well  as  those  published  in 
reference  9  are  plotted  in  Fig.  (8).  The  solid 
data  points  in  this  figure  represent  the  data  ob¬ 
tained  in  this  investigation  and  the  hollow  data 
points  represent  those  obtained  in  this  reference. 
From  this  figure  it  can  be  seen  that  the  relative 
changes  in  the  temperature  dependence  of  these 
three  aluminum  alloys  as  a  function  of  compress** 
stress  can  be  represented  by  Eq.  (2)  with  K  equa 
to  2.38  x  10"3  per  MPa.  This  result  emphasizes 
an  earlier  finding  concerning  the  insensitivity 
the  temperature  dependence  of  ultrasonic  velojjr, 
to  texture  and  alloy  composition,  and  makes  th 
method  more  suitable  for  the  nondestructive  «v* 
tion  of  stress  than  direct  velocity  measuremen  • 
The  latter  method  has  been  shown  to  be  strongiy^^ 
influenced  by  metallurgical  variables  and  req 
calibration  for  each  alloy  or  perhaps  for  each 
specimen1®. 

On  the  other  hand,  the  effect  of  stress  oj*^ 
temperature  dependence  of  ultrasonic  longitu 
velocity  becomes  much  smaller  and  opposite  * 
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.,en  the  stress  is  applied  parallel  to  theidiijk- 
,ion  i"  which  the  waves  are  propagated.  ThF^ 
alIrees  with  the  results  obtained  by  Chern  and 
u» yuan11  wb°  als0  observed  small  changes  in  the 
temperature  dependence  with  stress  in  2024 -T4 
aluminum  when  the  stress  was  applied  in  the  same 
direction  of  wave  propagation.  From  Figs.  6  and 

1  it  is  interesting  to  note  that  the  behavior  of 
the  temperature  dependence  with  stress  (Fig.  6) 
is  opposite  to  that  of  the  velocity  itself  (Fig. 

71.  The  comparison  between  the  effects  of  stress 
0n*the  temperature  dependence  and  on  the  velocity 
uill  be  discussed  below. 

2  Stress  Dependence  of  Ultrasonic  Velocity 

For  a  longitudinal  ultrasonic  wave  propagating 
parallel  to  the  direction  of  applied  stress,  the 
change  in  velocity  with  stress,  dV/do,  may  be  ex¬ 
pressed  as12 


where  f  is  the  round-trip  frequency  of  a  longitu¬ 
dinal  wave  propagating  in  the  axial  direction  of 
the  fastener,  V  is  the  ultrasonic  velocity,  E  is 
Young's  modulus  and  o  is  applied  stress.  For 
small  strains,  the  term  in  brackets  is  constant  and 
Eq.  (4)  can  be  written  as 


A  f  U-1 
-y  *  H  *  Ao 


dV  r  -1 

3^  2pV(3A  ♦  2u] 


[21*  X*  (4m ♦  4A  ♦ 

(3) 


In  Eq.  3  X  and  u  are  the  second  order  elastic  con¬ 
stants,  l  and  m  are  the  third  order  elastic  con¬ 
stants  of  Mumaghaml2t  p  js  density,  V  is  ultra¬ 
sonic  velocity,  and  o  is  applied  stress.  This 
equation  indicates  that  for  small  stresses,  dV/do 
is  a  constant.  Heyman  and  Chern have  used  the 
following  relationship  to  measure  axial  stresses 
in  fasteners 

6f  ,1  dv  ,,, 

T  lv  T'  r,Ao  (4) 


u-i  -  r  1  dV  1,-1 

H  -  l¥  TT  -  T1 


For  aluminum,  they  found  H  to  be  a  constant  equal 
to  -1.86  x  104  MPa. 

In  the  present  investigation,  the  effects  of 
compressive  applied  stress  on  ultrasonic  longitu¬ 
dinal  velocity  are  shown  in  Fig.  7.  From  this 
figure,  it  can  be  seen  that  for  stress  applied  in 
a  direction  parallel  to  that  of  ultrasonic  propa¬ 
gation,  the  results  give  a  value  of  -0.210  m/s-MPa 
for  the  stress  dependence,  dV/do.  Substituting  the 
values  V  -  6200  m/s,  E  -  70.3  GPa,  and  dV/do  =  -0.210 
m/s -MPa  into  Eq.  (6)  yields  H*1  «  -1.96  x  104  MPa. 
This  value  is  in  excellent  agreement  with  the  value 
of  H~1  equal  to  -1.86  x  104  MPa  obtained  by  Heyman 
and  Chern  for  2024-T6  aluminum. 

For  stress  applied  perpendicular  to  wave  pro¬ 
pagation.  the  change  in  velocity  with  stress  is 
given  by*2 


do  °  ~SpV(3A  ♦  2u)  l2t  '  ‘T  C"*^2*)] 

Table  (3)  -  Sunmiary  of  important  quantities  found  in  this  investigation  for  aluminum  6061-T6. 


Quantity 

‘#1  !  0  "  0 


[if]o  -  i-Sh 


.  3T  .  .  Room 
l  3o  -*T  ’  Temperature 


Stress  J.  Propagation 


Stress  //  Propagation 


-23.4  x 10 


0.0832 


4  »/s 


8.50  x  10 


-4  m/s 
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NONDESTRUCTIVE  STIXSS  MEASUREMENTS  IN  ALUMINUM 


K.  Salaaa,  G.  C.  Barber  and  N.  Chandrasekaran 
Department  of  Mechanical  Engineering 
University  of  Houston 
Houston,  Texas  77004 


Abstract 

The  effects  of  applied  elastic  stress  on  the  temperature  dependence  of 
10  Mil  ultrasonic  longitudinal  velocity  have  been  studied  in  aluminum  6061- 
T6.  Velocities  of  longitudinal  ultrasonic  waves  were  measured  as  a  function 
of  temperature  ir.  several  specimens  using  a  pulse -echo-overlap  system.  Ex¬ 
periments  were  performed  with  the  stress  applied  in  a  direction  parallel 
tc  and  perpendicular  to  the  ultrasonic  propagation  direction.  In  all 
temperature  dependence  measurements,  the  ultrasonic  velocity  is  found  to 
decrease  linearly  with  temperature,  and  the  slope  of  the  line  of  best  fit 
of  ultrasonic  velocity  versus  temperature  is  found  to  vary  considerably 
when  the  specimens  are  subjected  to  stress.  The  results  obtained  when  the 
stress  is  applied  in  a  direction  parallel  to  the  ultrasonic  propagation 
show  that  the  temperature  dependence  of  ultrasonic  velocity  increases 
linearl)  with  either  applied  tensile  or  compressive  stress.  In  the  case 
of  stress  applied  perpendicular  to  the  ultrasonic  propagation,  the  results 
indicate  that  the  te^>erature  dependence  decreases  linearly  with  either 
applied  tensile  ot  compressive  stress.  Calibration  curves  relating  the 
relative  change  in  the  temperature  dependence  of  ultrasonic  velocity  to 
applied  stress  are  constructed.  Using  these  calibration  curves,  the 
sensitivity  in  determining  unknown  applied  stresses  in  aluminum  is  estimat¬ 
ed  to  be  i  10  MPa. 


1 .  INTRODUCTION 


Materials  in  machine  components  are  always  in  a 
state  of  stress.  This  stress  can  be  applied 
stress  due  to  external  loading,  residual  stress 
within  the  material  without  any  external  load¬ 
ing,  or  a  combination  of  applied  and  residua] 
Applied  stresses  car.  be  calculated  using  formu¬ 
lae  from  mechanics  of  materials  or  can  be  found 
experimentally  using  strain  measuring  devices 
attached  to  the  machine.  Residual  stresses, 
however,  can  seldom  be  calculated  because  there 
is  usually  no  data  or.  which  stress  calculations 
are  based.  There  are  presently  several  methods 
of  experimentally  finding  residual  stresses  by 
destructive  means,  such  as  hole  drilling  and 
ring  coring2.  To  measure  residual  stresses  in 
a  body  non-destructivelv ,  however,  several  ap¬ 
proaches  based  on  three  major  methods  have  been 
proposed.  These  methods  aTe  x-ray  diffraction, 
electromagnetics,  and  ultrasonics.  None  of 
these  methods  measures  stress  directly.  X-ray 
diffraction  measures  lattice  strain,  electro¬ 
magnetic  methods  measure  magnetic  properties, 
and  ultrasonic  methods  aeasure  ultrasonic 


velocity.  All  of  these  methods  have  limitations 
which  prevent  any  of  them  from  being  used  in  all 
stress  measurement  applications.  For  the  non¬ 
destructive  evaluation  of  bulk  residual  stresses 
in  crystalline  and  non-crystalline  materials, 
the  methods  using  ultrasonic  techniques  seem  to 
hold  the  best  promise. 

The  temperature  dependences  of  the  elastic  con¬ 
stants  of  a  solid  are  due  to  the  enharmonic 
nature  of  the  crystal  lattice. A  measure  of 
the  temperature  dependence  of  the  ultrasonic  ve¬ 
locity  can,  therefore,  be  used  to  evaluate  bulk 
stresses.  Experiments  undertaken  on  aluminum 
and  copperS.6  elastically  deformed  in  compres¬ 
sion  showed  that  the  ultrasonic  velocity,  in 
the  vicinity  of  room  temperature,  changed  lin¬ 
early  with  temperature,  and  the  slope  of  the 
lineaT  relationship  changed  considerably  as  the 
amount  of  applied  stress  was  varied.  In  alumi¬ 
num,  the  relative  changes  of  the  temperature  de¬ 
pendence  of  longitudinal  velocity  decreased 
linearly  by  ss  much  as  23%  at  a  stress  of 


approximately  96  MPa.  The  lesultf  obtained  or. 
different  types  of  aluminum  alloys  also  indicate 
that  the  relative  changes  in  the  temperature 
dependence  due  to  applied  stress  are  insensi¬ 
tive  to  composition  and  texture,  and  the  data 
obatained  on  these  alloys  can  be  represented  by 
a  single  relationship7,  as 


where  (dV/dT)  is  the  temperature  dependence  at 
zero  applied  stress,  (dV/dT).  is  the  tempera¬ 
ture  dependence  at  an  applied  stress  c,  and  K 
is  a  constant  equal  to  2.3  x  10'3  per  MPa.  In 
addition,  the  effects  of  tensile  elastic  stress¬ 
es  on  the  temperature  dependence  of  longitudinal 
velocity  were  studied  using  the  aluminum  alloys 
2024-0,  3003-T251 ,  and  2024-T351  by  Salama  and 
Kang  (5) .  The  results  from  these  experiments 
were  found  to  satisfy  Eq.  (1)  with  the  constant 
K  equal  to  -1.89  x  10"3  per  MPa.  All  the  above 
studies  have  been'  performed  when  the  stress  was 
applied  in  a  direction  perpendicular  to  that 
of  the  ultrasonic  wave  propagation. 

In  the  present  investigation,  experiments  were 
performed  to  study  the  effects  of  both  compres¬ 
sive  and  tensile  elastic  stresses  on  the  tem¬ 
perature  dependence  of  longitudinal  ultrasonic 
velocity  in  three  specimens  of  the  aluminum  al¬ 
loy  6061-T6.  The  experiments  were  performed 
with  the  stress  applied  in  a  direction  parallel 
to  and  perpendicular  to  the  ultrasonic  wave  pro- 
pagal ion. 

2.  EXPERIMENTAL  PROCEDURE 


are  usee  in  the  present  investigation.  These 
arrangements  are  described  in  detail  else¬ 
where3  »fc.  To  help  ensure  the  uniformity  of 
stress  in  specimens,  special  effort  was  made  in 
designing  these  stress  application  systems  to 
minimize  any  effects  of  misalignment  between  the 
axis  of  the  specimen  and  the  loading  axis.  Por 
tensile  testing,  this  was  achieved  by  using  two 
universal  joints  between  specimen  and  coinec- 
tions  to  load  cell  and  to  the  loading  frame. 

For  compressive  testing,  a  two-stepped  alignment 
system  was  used.  A  linear  ball  bearing  served 
as  the  first  alignment  between  the  upper  and 
lower  loading  axes  and  the  hemispheric  steel 
ball  served  as  the  second  alignment  between  the 
specimen  and  the  loading  axis. 

The  ultrasonic  velocity  was  measured  using  the 
pulse-echo-overlap  method  which  has  been  fully 
described  elsewhere®.  X-cut  transducers  are 
used  for  the  generation  of  the  longitudinal 
waves,  which  were  used  in  all  measurements.  The 
transducer  is  placed  on  the  specimen  by  means  of 
special  holders  which  are  designed  to  clamp  to 
the  specimen  as  shown  in  Fig.  2  for  the  parallel 
configuration  and  Fig.  3  for  the  perpendicular 
configuration.  The  spring-supported  plunger 
serves  as  the  inner  conductor  of  the  coaxial 
cable  which  carries  the  10  MHz  electrical  sig¬ 
nals  to  and  fro r  the  transducer  which  is  bonded 
to  the  specimen.  These  signals  are  transmitted 
fror  the  plunger  by  a  piece  of  teflon-coated 
wire  to  a  BXC  receptacle  which  is  mounted  direct¬ 
ly  to  the  transducer  holder.  The  spring  is  used 
to  produce  a  pressure-coupled  transducer  which 
is  required  to  insure  uniformity  of  the  thin 
layer  bond.  A  small  clamping  force  of  30  to  SO 
Newtons  was  used  in  the  present  work. 


The  specimens  used  in  the  present  study  were 
made  from  one  inch  diameter  bar  stock  of  com¬ 
mercial  6061 -T6  aluminum  in  the  form  shown  in 
Fig.  1.  All  specimens  were  made  identical  ex¬ 
cept  the  overall  length  L  was  varied.  The 
specimens  were  machined  with  a  2.54  cm  diameter 
cap  on  each  end  which  allowed  the  same  speci¬ 
men  to  be  used  for  stress  applied  in  compres¬ 
sion  as  well  as  in  tension.  The  two  caps  at 
the  ends  were  made  flat  and  parallel  to  within 
i0.00:  cm  in  order  to  avoid  diffraction  effects 
in  the  ultrasonic  beam  during  propagation. 

These  caps  were  also  connected  to  the  center 
portion  of  the  specimen  by  a  0.06  cm  radius  in 
order  to  minimize  stress  concentrations.  AfteT 
experiments  of  the  stress  applied  parallel  to 
the  axis  of  specimen  were  completed,  two  paral- 
el  surfaces  of  thickness,  t,  were  milled  in  the 
center  of  the  specimen  to  allow  measurements 
for  stress  applied  perpendicular  to  wave  pro¬ 
pagation  direction.  In  this  case  the  ultra¬ 
sonic  waves  were  propagated  in  the  center  por¬ 
tion  of  the  specimen  where  the  stress  is  ex¬ 
pected  to  be  uniform. 


A  temperature  control  system  is  designed  to  en¬ 
close  the  specimen  gripping  assemblies  to  en¬ 
sure  stabilized  temperature  foT  the  entire  speci¬ 
men  during  the  time  required  for  velocity  measure 
ments.  All  experiments  were  performed  between 
approximately  310*li  and  370®K  where  the  coupling 
condition  between  the  transducer  and  the  speci¬ 
men  was  found  to  be  satisfactory.  A  test  fur¬ 
nace  was  used  for  heating.  This  furnace  pro¬ 
vides  a  uniform  temperature  performance  by  making 
use  of  a  recirculating  blower  and  planum  system. 
The  furnace  is  equipped  with  a  temperature  con¬ 
troller  which  is  capable  of  providing  precision 
temperature  control.  Before  taking  velocity 
measurements,  the  furnace  is  set  at  a  desired 
maximum  temperature.  After  this  desired  tem¬ 
perature  is  reached,  the  furnace  is  turned  off 
and  the  cooling  rate  of  the  specimen  is  con¬ 
trolled  to  achieve  a  constant  cooling  rate.  The 
temperature  of  the  specimen  is  measured  by  a 
copper- const ant an  thermocouple  attached  directly 
to  the  specimen.  The  thermocouple  along  with 
a  potentiometer  provides  the  measurement  of  the 
specimen  temperature  to  an  accuracy  tO.l'K. 


The  application  of  external  stress  was  carried  3.  RESULTS 

out  with  a  model  1125  floor  type  lnstron 
machine  of  20,000  kg  maximum  load  capacity. 

Four  different  types  of  loading  arrangements 


The  velocity  of  longitudinal  ultrasonic  waves 
propagating  perpendicular  to  the  applied  stress 


hat  r.casurtc  at  a  function  of  temperature  in 
specimens  3  and  2.  Typical  examples  of  the  ul¬ 
trasonic  velocity  obtained  on  specimen  2  at  the 
applied  tensile  C.O  MPa  and  66.6  MPa  are  shown 
in  Fig.  4.  This  figure  shows  that,  in  the  tem¬ 
perature  range  from  330*K  to  370'K,  ultrasonic 
velocity  decreases  linearly  with  temperature. 

A  computer  program  was  used  to  determine  the 
temperature  dependence  of  ultrasonic  velocity, 
dV/dT. 


for  the  experimental  results  are  given  in  Table 
(2).  Also  listed  in  Table  (1)  are  the  relative 
changes  in  the  temperature  dependence,  due  to 
application  of  stress.  These  relative  percentage 
changes  of  the  temperature  dependence  are  plot¬ 
ted  as  a  function  of  applied  stress  in  Fig. 

(S).  From  this  figure,  it  can  be  seen  that  the 
data  points  for  tensile  stress  lie  on  a  straight 
line  which  passes  through  the  origin.  This  line 
has  a  slope  of  -1.16  x  10'^ /MPa  and  a  correla¬ 
tion  coefficient  equal  to  -0.972.  The  data 


Table  (1)  -  Variations  of  the  temperature  dependence  of  longitudinal 
ultrasonic  velocity  with  applied  tensile  and  compres¬ 
sive  stress  in  6061-T6  aluminum.  The  stress  is  applied 
perpendicular  to  the  direction  of  ultrasonic  propagation. 


Specimen 

Number 

Specimen 
Dimensions 
o  finch) 

Applied 

Stress 

(MPa) 

dV/dT 

(m/s-K) 

Offset 

Uncertainty 

(m/s) 

L\ 

1 

t 

•  0.624  (0.24SS) 

66.2 

-.9177 

.158 

-8.04 

L 

«  6.985  (2.750) 

44.1 

-.9516 

.198 

-4.62 

22.3 

-.9530 

.304 

-4. SO 

0.0 

-1.0070 

.358 

0.00 

-16.6 

-.9657 

.326 

-3.86 

-33.1 

-.9333 

.358 

-7.11 

-66.1 

-.8631 

.197 

-14.09 

2 

t 

«  0.495  (0.19S) 

66.6 

-.9396 

.224 

-8.44 

L 

»  6.351  (2.5005) 

46.7 

-.9685 

.234 

-5.62 

26.7 

-.9973 

.195 

-2.81 

0.0 

-1.0230 

.192 

0.00 

Table  (2)  -  Values  of  slope,  Y-intercept,  and  correlation  coefficient  in  the  lines 
of  best  fit  of  dV/dT  versus  applied  stress  in  aluminum  6061 -T6.  The 
stress  is  applied  perpendicular  to  the  direction  of  ultrasonic  propa¬ 
gation. 


Specimen 

Number 

Specimen 
Dimensions 
cm  (inch) 

Applied 

Stress 

Slope 

(m/s-K  MPa) 

Correlation 

Coefficient 

Y-intercept 

(m/s-K) 

1 

t  «  0.624  (0.2455) 

L  *  6.985  (2.750 

Tensile 

1.22  x  JO'1* 

.941 

-0.9979 

2 

t  «  0.495  (0.195) 

L«  6.351  (2.5005) 

Tensile 

1.26 x 30’3 

.995 

-3.0262 

1 

t«  0.624  (0.2455) 
L«  6.985  (2.750) 

Compressive 

-2.15 x 10'3 

-.999 

-3.0047 

Table  (1)  lists  the  values  of  the  temperature 
dependence  of  longitudinal  velocity  and  the  of¬ 
fset  uncertainty  in  specimens  1  and  2  when  they 
were  subjected  to  applied  stresses  ranging  from 
•66.2  MPa  to  66.8  MPa.  This  table  shows  that 
the  magnitude  of  the  temperature  dependence  of 
ultrasonic  velocity  decreases  linearly  with  ap¬ 
plied  compressive  or  tensile  stress  with  the 
maximum  value  for  the  temperature  dependence  of 
ultrasonic  velocity  occuring  at  0.0  MPa.  The 
values  of  the  slope,  the  intercept,  and  the  cor¬ 
relation  coefficient  of  the  lines  of  best  fit 


points  for  compressive  stress  also  lie  on  a 
straight  line  which  passes  through  the  origin 
and  has  a  slope  of  2.44  x  10‘^/MPa  and  a  cor¬ 
relation  coefficient  equal  to  0.982. 

The  velocity  of  longitudinal  ultrasonic  waves 
propagating  parallel  to  applied  stress  was 
measured  as  a  funtion  of  temperature  on  the 
two  aluminum  specimens  2  and  3.  Table  (3) 
lists  the  results  obtained  on  these  two  alumi¬ 
num  6061 -T6  specimens  investigated  in  the 
parallel  configuration.  Included  in  this  table 


»Tt  the  value*  of  tni  temperature  dependence  of  tenperature  dependence  of  ultrasonic  velocity 

longitudinal  ultrasonic  velocity  determined  at  increases  linearly  with  either  applied  compres- 

applied  compressive  and  tensile  stresses  ranging  sive  or  tensile  stress,  with  the  minimum  value 

between  -79.5  Mia  and  S9.6  MPa.  Also  listed  in  of  the  temperature  dependence  occurring  at  0.0 

table  (3)  are  the  offset  uncertainities  which  MPa.  The  values  of  the  slope,  the  intercept, 

measure  the  linearity  of  the  relationship  and  the  correlation  coefficient  for  the  lines  of 

between  velocity  and  temperature.  The  results  best  fit  for  the  experimental  results  shown  in 

of  Table  (3)  are  plotted  in  Fig.  (6)  which  il-  Fig.  (6)  are  given  in  Table  (4).  From  Tables 

lustrates  the  effect  of  applied  co^>ressive  and  (2)  and  (4),  it  can  be  seen  that  the  slopes  of 

Table  (31  -  Variations  of  the  temperature  dependence  of  longitudinal  ultrasonic 
velocity  with  applied  tensile  and  compressive  stress  in  aluminum 
6061 -T6.  Stress  is  applied  parallel  to  the  direction  of  ultra¬ 
sonic  propagation. 


Specimen 

Number 

Specimen  Length 
cm  (inch) 

Applied  Stress 
(MPa) 

dV/dT 

(m/s-K) 

Offset 

Uncertainty 

(m/s) 

t\ 

2 

6.351  (1.5005) 

39.7 

-1.0529 

.423 

1.37 

29.8 

-1.0484 

.185 

.93 

19.8 

-1.0469 

.163 

.79 

9.9 

-1.0417 

.248 

.29 

0 

-1.0146 

.243 

0 

-19.8 

-1.0231 

.236 

.88 

-49.7 

-1.0355 

.266 

2.10 

-59.6 

-1.0427 

.351 

2.81 

3 

3.807  (1.499) 

59.6 

-1.0071 

.466 

3.14 

49. 7 

-.9960 

.318 

2.01 

29.8 

- .9944 

.530 

1.84 

9.9 

-.9798 

.326 

.35 

0 

-.9692 

.166 

0 

-19.8 

-.9820 

.234 

1.53 

-39.7 

-.9961 

.252 

2.99 

-59.6 

-1.0158 

.233 

5.02 

-79.5 

-1.0317 

.342 

6.69 

Table  (4)  -  Values  of  slope,  intercept,  and  correlation  coefficient  taken  from 

lines  of  best  fit  of  dV/dT  versus  applied  stress  in  aluminum  6061-T6. 
The  stress  is  applied  parallel  to  the  direction  of  ultrasonic  propa- 
gat ion . 


Spec imen 
Number 

Specimen 

Length 
cm  (inch) 

Applied 

Stress 

Slope 

(m/s*K  MPa) 

Correlation 

Coefficient 

Y- Intercept 
(m/s-K) 

1 

6.351  (2.5005) 

Tensile 

-3.53  x  lO-4 

-.981 

-1.0387 

3 

3.807  (1.499) 

Tensile 

-4.82 x 10"4 

-.948 

-.9764 

1 

6.351  (2. 5005) 

Compressive 

4.56 x 10"4 

996 

-1.0142 

3 

3.807  (1.499) 

Compressive 

7.99 x 10'4 

.997 

-.9672 

tensile  stresses  on  the  temperature  dependence 
of  ultrasonic  velocity.  This  figure  shows 
that,  within  the  applied  stress  range  used  in 
these  measurements,  the  magnitude  of  the 


the  lines  of  best  fit  for  the  perpendicular 
configuration  are  much  larger  and  opposite 
in  sign  to  those  obtained  in  the  parallel  con¬ 
figuration. 


Fjfure’*t  indicate*  that  the  plot  of  dV/d7 
versus  stress  for  a  gntr.  specimen  has  an  abrupt 
change  ir.  direction  where  the  applied  stress 
equals  zero.  The  intercept  foT  the  tensile 
portion  of  the  graph  is  also  significantly  high¬ 
er  than  the  intercept  for  the  compression  por¬ 
tion  of  the  graph.  In  the  case  of  tensile 
loading,  bending  stresses  are  introduced  and 
likely  to  exist  throughout  the  entire  length 
of  the  specimen.  This  non-uniform  loading  is 
expected  to  affect  the  absolute  magnitude  of 
dV/dT,  but  not  the  relative  differences  between 
quantities  as  indicated  by  the  consistency  of 
the  slopes  of  the  lines  of  best  fit  of  dV/dT 
versus  stress  in  these  specimens.  In  the  case 
of  compressive  loading,  the  bottom  end  cap  of 
the  specimen  is  loaded  uniformly  which  gives 
the  compression  loading  a  better  probability 
of  uniform  stress  throughout  the  ultrasonic 
path. 

Because  the  values  of  the  intercept  of  the  lines 
of  best  fit  were  found  to  vary  among  the  speci¬ 
mens  investigated,  the  relative  change  in  the 
temperature  dependence,  £,  due  to  the  applica¬ 
tion  of  stress  was  calculated  and  its  values  are 
listed  in  column  6  of  Table  (3).  The  values  of 
L  were  calculated  using  the  relationship, 

(dV/dT)  -  (dV/dT) 

L  =  (3v73t]  (2) 

where  (d\ /dT)  is  the  temperature  dependence 
given  by  the  intercept  of  the  line  of  best  fit 
and  (dV/dT}c  is  the  temperature  dependence  at 
ar.  applied  stress,  c.  The  variations  in  the 
temperature  dependence  of  the  three  specimens 
tested  at  zero  stress  are  believed  to  be  due  to 
differences  in  the  residual  stress  in  these 
specimens. 

The  relative  changes  in  the  temperature  depen¬ 
dence,  obtained  or.  all  three  specimens  investi¬ 
gated,  are  plotted  in  Fig.  7.  The  plot  for 
tensile  stress  shows  that  the  points  lie  on  a 
straight  line  which  passes  through  the  origin 
with  a  slope  equal  to  4.39 x  10~4^agr  MPa  and  a 
correlation  coefficient  equal  to  0.930.  The 
plots  for  compressive  stress  show  that  the 
points  for  specimen  3  lie  on  a  straight  line 
that  passes  through  the  origin  with  a  slope 
equal  to  -8.63  x  10*4  per  MPa  and  a  correlation 
coefficient  equal  to  -0.992.  The  points  for 
specimen  2  for  compressive  stress  also  lie  on 
a  straight  line  that  passes  through  the  origin 
with  a  slope  equal  to  -4.56  x  10"4  per  MPa  and  a 
correlation  coefficient  equal  to  -0.996.  There 
is  no  explanation  for  why  the  slope  of  the  data 
obtained  on  specimen  2  is  smaller  than  that 
found  from  the  data  obtained  on  specimens  3. 

4.  DISCUSSION 

In  the  present  investigation,  experiments  were 
performed  to  study  the  effects  of  both  compres¬ 
sive  and  tensile  elastic  stresses  on  the  temp¬ 
erature  dependence  of  longitudinal  ultrasonic 


velocity  measured  in  a  direction  perpendicular 
to  that  of  stress  in  three  specimens  of  the 
aluminum  alloy  6061-T6.  The  relative  changes 
in  the  temperature  dependence  of  longitudinal 
ultrasonic  velocity  as  a  function  applied  stress 
obtained  in  this  work  as  well  as  those  published 
in  references  7  and  9  are  plotted  in  Fig.  (8). 

The  solid  data  points  in  this  figure  represents 
the  results  obtained  in  this  investigation  and 
the  hollow  data  points  represent  those  obtained 
in  these  references.  From  this  figure  it  can 
be  seen  that  the  relative  changes  in  the  tempera¬ 
ture  dependence  of  alumima  alloys  as  a  function  of 
compressive  stress  can  be  represented  by  Eq.  (1) 
with  K  equal  to  2.38 x 10~3  per  MPa.  The  line  of 
best  fit  of  the  relative  changes  in  temperature 
dependence  versus  tensile  stress  for  aluminum 
6061-T6  is  seen  to  vary  slightly  from  the  line 
of  best  fit  of  the  other  aluminum  alloys.  This 
variance  is  believed  to  be  due  to  the  non-uniform 
stress  induced  in  the  specimens  in  the  case  of 
tensile  loading.  Therefore,  the  line  of  best 
fit  for  tensile  stress  for  the  aluminum  alloys 
2024-0,  3003-T2S1,  and  2024-T351  is  chosen  to 
represent  the  true  line  of  best  fit  of  aluminuir 
alloys. 

In  OTder  to  use  the  temperature  dependence  method 
to  determine  unknown  applied  stresses,  a  calibra¬ 
tion  curve  was  constructed  and  is  shown  in  Fig. 

9.  This  calibration  curve  shows  the  lines  of 
best  fit  of  the  relative  percentage  changes  of 
the  temperature  dependence  versus  applied  stress 
found  in  Fig.  (8).  The  bands  drawn  around  these 
lines  of  best  fit  are  for  a  95»  confidence  level 
for  stress.  These  bands  were  determined  by  means 
Of  an  inverse  regression  technique  which  is  des¬ 
cribed  in  detail  in  reference  10.  From  Figure 
(9)  it  can  be  seen  that  the  relative  change  in 
temperature  dependence  decreases  when  either 
compressive  or  tensile  stresses  are  applied  to 
the  specimen.  No  explanation  is  available  at 
present  for  this  behavior  which  makes  it  diffi¬ 
cult  to  differentiate  between  tensile  and  compres¬ 
sive  stresses  when  using  the  temperature  depen¬ 
dence  method  in  the  non-destructive  evaluation 
of  applied  stress  in  aluminum  alloys.  To  use 
the  calibration  curve  shown  in  Fig.  (9)  to  deter¬ 
mine  applied  stresses,  the  temperature  dependence 
of  ultrasonic  velocity  is  measured  in  the  speci¬ 
men  at  zero  applied  stress.  This  measurement 
should  be  repeated  several  times  and  the  average 
value  is  taken  as  the  true  value.  The  specimen 
is  then  loaded  and  the  temperature  dependence 
is  measured  perpendicular  to  the  applied  stress. 
Using  the  values  found  above,  the  relative  per¬ 
centage  change  in  the  temperature  dependence  of 
ultrasonic  velocity  if  calculated  using  Eq.  1- 
This  value  along  with  Fig.  (9)  are  used  to  deter¬ 
mine  the  applied  stress.  If  greater  accuracy  is 
required,  the  temperature  dependence  at  the  un¬ 
known  stress  should  be  measured  several  times  and 
the  average  value  used  to  determine  the  applied 
stress. 


The  effects  of  compressive  and  tensile  stress  on 
the  temperature  dependence  of  longitudinal  ul¬ 
trasonic  velocity  were  measured  on  two  specimens 
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(2  and  3;  os  aluminur  OObl-Te.  These  experiment; 
weTt  pej forced  with  the  stress  applied  ir  a 
direction  parallel  to  the  ultrasonic  propaga¬ 
tion.  These  experiments  showed  that  the  relative 
change  in  the  temperature  dependence  of  ultrasonic 
velocity  is  a  linear  function  of  the  applied 
stress  and  can  be  represented  by  Eq.  (1).  The 
values  of  K  found  from  fig.  (7)  are  4.39xl0*4 
per  MPa  for  tensile  stresses,  -4.56xl0"4  per 
MPa  for  specimens  2  in  compression,  and 
-8.63  x  10~4  per  MPa  for  specimen  3  in  compres¬ 
sion. 

In  order  to  use  the  temperature  dependence 
method  to  determine  unknown  applied  stress,  a 
calibration  curve  was  constructed  and  is  shown 
in  Fig.  10.  From  this  figure  it  can  be  seen 
that  the  relative  change  in  temperature  depen¬ 
dence  increases  when  either  compressive  or 
tensile  stress  is  applied  to  the  specimen. 

This  makes  it  difficult  to  differentiate  be¬ 
tween  tensile  and  compressive  stress  (as  was 
the  case  in  the  perpendicular  configuration) 
when  using  the  temperature  dependence  method  in 
the  non-destructive  evaluation  of  applied  stress 
in  aluminum.  6061-T6. 

To  demonstrate  the  use  of  this  calibration  curve 
consider  an  example  where  the  temperature  de¬ 
pendence  is  measured  parallel  to  an  applied 
compressive  stress  and  the  relative  percentage 
change  in  the  temperature  dependence  calculated 
to  be  3.04.  The  accuracy  of  dV/dT  measurements 
in  these  experiments  was  found  to  be  11.64 
using  the  method  of  standard  error  estimation.11 
Thus  in  this  example,  the  true  value  of  the 
temperature  dependence  will  be  311.64.  The 
applied  stress  is  read  from  Fig.  do)  as  -35 
MPa  with  a  lower  bound  of  -57  MPa  and  an  uppeT 
bound  of  -13  MPa.  This  amount  of  variance  is 
large  and  would  be  unacceptable  for  many  appli¬ 
cations.  This  variance,  howeveT.  can  be  re¬ 
duced  by  measuring  the  temperature  dependence 
at  the  applied  compressive  stress  several  times. 
Consider  the  case  where  the  temperature  depen¬ 
dence  at  an  applied  compressive  stress  is 
measured  five  times  and  the  average  value  is 
used  to  calculate  a  relative  percentage  change 
of  (dV/dT)  of  3.004.  The  standard  error  in 
this  example  is  1.6//5  *  0.74.  The  true  value 
of  the  relative  percentage  of  (dV/dT)  is 
3.00  1  .74.  The  applied  stress  is  read  from 
Fig.  10  as  -35  MPa  with  a  lower  bound  of  -24 
MPa  and  an  upper  bound  of  -46  MPa.  These  ex¬ 
amples  show  that  the  temperature  dependence 
measurement  should  be  repeated  several  times 
when  the  stress  is  applied  in  a  direction 
parallel  to  the  ultrasonic  propagation.  This 
will  produce  a  stress  measurement  with  an 
acceptable  amount  of  error  for  most  applica¬ 
tions. 
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Transducer  holder  used  for  the  measure¬ 
ment  of  ultrasonic  velocity  when  the 
stress  is  applied  in  a  direction  perpen¬ 
dicular  to  ultrasonic  propagation. 


Fig.  2  Transducer  holder  used  for  the  measure¬ 
ment  of  ultrasonic  velocity  when  the 
stress  is  applied  in  a  direction  parallel 
to  ultrasonic  propagation. 
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Fig.  4  Effect  of  applied  tensile  stress  on  the 
temperature  dependence  of  ultrasonic 
longitudinal  velocity  in  alroinum  alloy 
6061 -T6. 


Relative  percentage  change  in  the 
temperature  dependence  of  ultrasonic 
longitudinal  velocity  as  a  function  of 
applied  compressive  and  tensile  stresses 
in  aluainw  6061-T6. 


Fig.  6  Effect  of  applied  compressive  and  tensile 
stresses  on  the  temperature  dependence 
of  ultrasonic  longitudinal  velocity. 
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Fig.  7  Relative  percentage  change  in  the 

tenperature  dependence  of  ultrasonic 
longitudinal  velocity  as  a  function  of 
applied  compressive  and  tensile  stress 
in  aluminum  6061-T6. 
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Fig.  E  Relative  percentage  change  in  the 

temperature  dependence  of  longitudinal 
ultrasonic  velocity  as  a  function  of 
applied  stress  in  aluminum  alloys.  The 
'hollow'  data  points  were  obtained  by 
Salama,  Ling  and  hang7 >9. 


Fig.  9  Calibration  curve  using  the  relative  Fig.  10  Calibration  curve  using  the  relative 

change  of  the  temperature  dependence  change  of  the  temperature  dependence 

of  ultrasonic  velocity  to  determine  of  ultrasonic  velocity  to  determine 

applied  stress  in  aluminum  alloys.  applied  stress  in  aluminim  6061-T6. 
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SYMPOSIUM  ON  APPLICATIONS  AND  DEVELOPMENT  OF  NONDESTRUCTIVE 
EVALUATION  FOR  USE  IN  MATERIALS  PROCESSING  -  OCT.  1983 

DETERMINATION  OF  STRESS  GENERATED  BY  SHRINK  FIT 

N.  Chandrasekaran,  H.  Y.  Wu  and  K.  Salama 
University  of  Houston 


The  temperature  dependence  of  ultrasonic  ve¬ 
locity  was  used  to  determine  the  stress  generat¬ 
ed  in  a  disc  made  of  AS33B  steel.  The  stress  is 
introduced  by  shrink  fitting  a  rod  of  the  same 
steel  in  a  hole  of  slightly  smaller  size.  The 
measurements  were  made  using  longitudinal  and 
shear  waves  propagating  along  the  thickness  of 
the  disc  in  the  temperature  range  between  room 
temperature  and  -40*C.  In  addition,  shear 
velocities  were  measured  with  the  polarization 
adjusted  parallel  to  and  perpendicular  to  the 
radius  on  which  the  stress  measurements  were 
performed.  The  results  show  that  at  all 
locations,  the  ultrasonic  velocity  decreases 
linearly  as  the  temperature  is  increased  and 
that  the  temperature  dependence  of  ultrasonic 
velocity  can  be  determined  to  J2%.  The  results 
also  indicate  that  the  temperature  dependence 
of  the  longitudinal  velocity  remains  unchanged 
as  the  radial  distance  on  the  disc  is  changed, 
which  agrees  with  predictions  of  elasticity 
theory.  The  temperature  dependences  of  the 
shear  velocities,  however,  show  considerable 
variations  with  radial  distance,  reflecting 
the  changes  in  the  radial  and  the  tangential 
components  of  stress  as  a  function  of  the 
radial  distance. 


ULTRASONIC  METHODS  SEEM  TO  HOLD  Tills  BEST 
promise  for  the  nondestructive  evaluation  of 
bulk  residual  stresses  in  both  crystalline 
and  non-crystalline  materials^.  Several 
methods  based  on  the  measurements  of  ultrasonic 
velocities?.3  or  their  temperature  dependences4 
are  currently  being  developed  and  show  promis¬ 
ing  possibilities  for  their  applications  in 
stress  measurements.  In  order  for  these 
methods  to  provide  quantitative  determination 
of  stress,  specimens  with  well-known  stress 
distributions  should  be  available  for  refer¬ 
ence.  Ilsu  et  alj>  have  studied-  the  distribu¬ 
tion  of  ultrasonic  longitudinal  velocity  in  a 
shrink-fit  assembly  and  found  a  good  agreement 


between  the  velocity  measurements  and  the  cal¬ 
culated  stress  distribution. 

The  temperature  dependences  of  ultrasonic 
velocities  are  due  to  the  anharmonic  nature  of 
the  solid,  and  therefore  can  be  used  for  stress 
measurements.  Experiments  performed  on  alumi¬ 
num6,  copper  and  steel7  show  that,  in  the 
vicinity  of  room  temperature,  the  ultrasonic 
longitudinal  velocity  decreases  linearly  with 
temperature,  and  the  slope  of  the  linear  re¬ 
lationship  varies  considerably  when  the  specimen 
is  subjected  to  stress.  The  results  also, 
indicate  that  the  temperature  dependence  of  ul¬ 
trasonic  velocity  is  insensitive  to  variations 
in  alloy -chemical  composition  which  represents 
a  major  difficulty  for  the  use  of  direct  velocity 
measurements  in  the  determination  of  residual 
stresses6. 

In  order  to  further  develop  the  temperature 
dependence  method  for  field  applications,  the 
radial  distribution  of  the  slope  of  ultrasonic 
velocity  vs.  temperature  has  been  studied  in  a 
Shrink-fit  assembly  made  of  ASTM  A533B  steel. 

The  study  is  made  using  longitudinal  and  shear 
waves  polarized  parallel  as  well  as  perpendicular 
to  radial  direction.  The  results  show  that  the 
temperature  dependence  Df  ultrasonic  velocity  re¬ 
mains  unchanged  throughout  the  radius  of  the 
stressed  specimen,  which  agrees  with  the  predic¬ 
tions  of  elasticity  calculations.  The  behavior 
of  the  temperature  dependences  of  shear  waves, 
however,  is  found  to  have  less  agreement. 

EXPERIMENTAL 

SPECIMEN  PREPARATION  -  Carbon  steel  ASTM 
A533B  was  used  in  the  present  research.  This 
steel  is  generally  used  in  the  manufacture  of 
pressure  vessels  in  the  nuclear  industry.  A 
specimen  of  ring-plug  type  was  machined  from 
this  steel.  The  dimensional  details  of  the  ring 
specimen  is  shown  in  Fig.  1.  The  figure  shows 
the  difference  in  the  dimension  of  outer  diameter 
of  the  plug  and  the  inner  diameter  of  the  ring 
before  plugging.  This  difference  is  the  extent 
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Fig.  1  Dimensional  details  of  reference 
specimen  »l. 


of  diametrical  interference.  The  mating  sur¬ 
faces  of  this  specimen  was  carefully  machined 
to  within  0.0005  in.  tolerance  which  was  met 
at  only  three  fourths  of  the  circumference. 

At  the  worst  point  the  deviation  from  cir¬ 
cularity  was  0.002  in. 

The  interference  fit  was  made  by  simul¬ 
taneously  cooling  the  plug  and  heating  the 
ring,  and  thereby  creating  a  temperature  dif¬ 
ferential  enough  to  Compensate  for  the  diff¬ 
erence  in  dimension  between  the  ring  and  the 
plug.  The  plug  was  immersed  in  liquid  nitrogen 
for  about  5  minutes  while  the  ring  was  heated  to 
about  60U*F  for  about  half  an  hour.  This 
temperature  difference  was  found  to  be  adequate 
for  plugging. 

TEMPERATURE  CONTROL  SYSTEM  -  The  temperature 
control  system  is  designed  to  enclose  the  speci¬ 
men  and  the  measuring  accessories  to  ensure  a 
stablized  temperature  during  the  time  required 
for  the  velocity  measurement.  The  desired 
subambient  temperature  is  obtained  by  control¬ 
ling  the  flow  rate  of  liquid  nitrogen  to  the 
insulated  enclosure.  After  the  specimen  is 
cooled  to  the  desired  low  temperature,  the 
supply  of  nitrogen  is  stopped.  The  warming  rate 
is  adjusted  so  that  the  rise  in  temperature  is 
steady  at  about  0.25*K/minute.  This  rate  is 
achieved  by  using  a  low  heat  input  light  source 
within  the  enclosure.  All  the  measurements  were 
made  in  the  range  of  220*K  to  260*K,  which  pro¬ 
duced  a  steady  warming  rate  for  the  system.  The 
actual  temperature  of  the  specimen  was  measured 
by  means  of  a  Copper-Constantan  thermocouple 
attached  directly  to  the  specimen.  The  thermo¬ 
couple  along  with  a  potentiometer  provides  the 
measurement  of  the  temperature  to  an  accuracy  of 
tO. IK . 


VELOCITY  MEASUREMENT  SYSTEM  -  Velocity 
measurements  in  the  present  investigation  was 
performed  by  the  pulse-ccho-ovcrlap  method, 
which  has  been  fully  described  by  Papadakis.  A 
block  diagram  of  the  apparatus  is  shown  in  Fig. 
2.  In  this  system,  a  pulse  of  approximately  1 
microsecond  in  duration  and  of  variable  pulse 
repetition  rate  is  generated  by  the  r.f.  pulsed 
oscillator  and  is  impressed  on  a  0.25  inch  in 
diameter  quartz  transducer  of  a  fundamental 
frequency  of  10  Mlz.  The  transducer  is  acousti¬ 
cally  bonded  to  the  specimen  by  means  of  a  thin 
layer  of  viscous  Nonaq  stopcock  grease. 


Fig.  2  Pulse-echo-overlap  system. 


After  a  conventional  display  of  echoes  is 
displayed  on  the  screen  of  the  oscilloscope, 
two  echoes  are  visually  superimposed  by  adjust¬ 
ing  the  frequency  of  the  CW  oscillator  so  that 
the  period  of  this  frequency  is  exactly  equal  to 
the  round  trip  time  of  the  ultrasonic  pulse  in 
the  specimen.  The  frequency  f  at  the  exact 
cycle  for  cycle  match,  is  then  recorded  directly 
from  the  electronic  counter  and  used  to  compute 
the  ultrasonic  velocity  using  the  relationship 
V  •  21f,  where  21  is  the  total  path  length. 

Figure  3  shows  the  transducer  holder  used 
in  the  measurements  of  velocity.  The  spring 
supported  plunger  served  as  the  inner  conductor 
of  the  coaxial  cable  which  carries  the  rf  sig¬ 
nals.  These  signals  are  transmitted  from  the 
plunger  by  a  piece  of  teflon  coated  wire  to  a 
BNC  receptacle  mounted  directly  on  the  trans¬ 
ducer  holder.  A  clamping  force  of  30  to  50 
Newtons  is  found  to  be  sufficient  to  produce 
a  uniform  thin  layer  bond  between  transducer 
and  specimen. 


The  velocities  of  longitudinal  and  shear 
waves  propagating  through  the  thickness  of  the 
specimen  are  measured  at  various  locations 
along  that  radial  reference  line.  The  measure¬ 
ments  of  the  shear  velocity  arc  made  when  the 
waves  are  polarized  parallel  as  well  as  per¬ 
pendicular  to  the  reference  line.  The  results 
of  these  measurements  are  representative  to  the 
stress  distribution  in  the  entire  shrink-fit 


Fig.  3  Transducer  holder  used  for  the 

measurement  of  ultrasonic  velocity  on 
reference  specimens. 


Fig.  4  Specimen  used  in  the  ultrasonic  velocity 
measurements. 


specimen  since  the  triaxial  stress  generated 
by  the  interference  fit  is  radially  syanetric. 
The  velocity  measurements  are  also  performed 
at  various  stages  of  manufacturing  of  the 
specimen.  Shear  waves  polarized  parallel  to 
circumferential  direction  is  hereinafter  called 
shear-parallel  and  shear  waves  polarized  along 
the  radial  direction  is  called  shear-perpen¬ 
dicular,  as  shown  in  Fig.  4.  The  measurements 
are  made  at  about  six  locations,  along  the 
reference  radial  line. 

The  results  of  typical  examples  of  the 
temperature  variations  of  longitudinal,  shear- 
parallel  and  shear-perpendicular  velocities  on 
specimen  with  the  hole  plugged  are  shown  in 
Table  (1).  These  measurements  are  made  at  a 
distance  of  0.2S  in.  from  the  edged  of the  hole. 
Table  (1)  also  shows  the  velocity  deviation 
from  the  best  fit  line  and  the  values  of  the 
dV/dT.  Fig.  S  shows  plots  of  the  results  of 
these  three  examples.  From  this  figure  it  can 
be  seen  that  the  velocity  of  ultrasonic  waves 
varies  linearly  with  temperature  and  the  slopes 
of  the  straight  line  of  best  fit  are  different. 

The  velocities  of  longitudinal  and  shear- 
parallel  and  shear-perpendicular  are  measured 
as  a  function  of  temperature  before  and  after 
stresses  are  introduced  in  specimen.  These 
measurements  are  made  along  the  sane  radial 
reference  line.  The  results  of  the  temperature 
dependences  of  ultrasonic  velocities  prior  to 
making  the  interference  fit  between  the  plug 
and  the  ring,  form  the  baseline  data  for  the 
study  of  stresses  generated  by  shrink  fit. 
Separate  baseline  data  is  obtained  for  the 
longitudinal,  the  shear-parallel,  and  the  shear- 
perpendicular  velocities  for  the  specimen. 

Figure  6  shows  the  variation  of  the  tem¬ 
perature  dependence  of  ultrasonic  longitudinal 
velocity  as  a  function  of  radial  distance  from 
the  hole.  The  plot  shows  the  variation  of 
(dV/dT) g  before  and  after  stresses  are  intro¬ 
duced.  The  results  of  the  temperature 


dependences  of  ultrasonic  shear  velocities  for 
the  same  two  cases  are  shown  as  a  function  of 
radial  distance  from  the  hole  in  Figs.  7  and  8. 
Figure  7  shows  the  variation  of  (dV/dT)j,  while 
Fig.  8  shows  the  variation  of  (dV/dT)^.  In  all 
cases  the  ultrasonic  waves  are  propagated  in 
the  axial  direction. 

DISCUSSION 

ACCURACY  OF  MEASUREMENTS  -  In  order  to 
estimate  the  repeatability  and  the  accuracy  of 
the  measuring  technique  the  measurements  of  ul¬ 
trasonic  velocity  as  a  function  of  temperature 
were  repeated  for  the  cases  of  longitudinal, 
Shear-parallel  and  shear-perpendicular.  The 
experiments  were  repeated  on  several  different 
dates  using  different  wave  propagation  nodes  and 
polarizations.  The  results  of  these  experiments 
are  tabulated  in  Table  (2).  The  table  also  shows 
the  standard  deviation,  p,  based  on  the  results 
of  the  four  sets  of  measurements  performed.  The 
standard  error,  S_,  is  computed  using  the  rela¬ 
tionship  of  eqn  (T),  to  obtain  the  error  in 
determining  the  temperature  dependence  of  ultra¬ 
sonic  velocity,  dV/dT® 

SB  ■  p//n  (1) 

where  n  is  the  number  of  measurements  and  p  is 
the  standard  deviation  of  an  infinite  number 
of  measurements. 

In  this  investigation,  dV/dT  is  measured 
only  once  for  each  particular  type  of  wave  pro¬ 
pagation.  Therefore,  equation  (1)  becomes. 


where  Sgg  is  the  standard  error  for  single 
measurement.  Using  the  results  in  table  (2) 
and  equation  (2),  S  t  is  found  to  be  £1.7%. 
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Table  (1) 

-  Variations 

of  ultrasonic 

longitudinal,  shear  (parallel)  and  shear  (perpendicular) 

velocities 

with  temperature  in  reference  specimen  #1,  after  the  interference  fit. 

/  •*.  1 

Measurements  were  done  at 

the  same  location 

of  0.25  in.  from  the  hole  edge  and  the 

wave  propagation  it  from 

surface  to  surface 

— ■  ^ — — 1 

ft  . 

Pol. 

Temperature 

Frequency 

Velocity 

Vel .  Dev. 

Deg  K 

Hz 

m/sec 

m/sec 

Long 

215 

232889 

5927.2 

-0.05 

«*  *-•  'j 

220 

232794 

5924.8 

0.28 

222 

232749 

5923.7 

0.11 

227 

232654 

5921.2 

-0.09 

ft  ,.-M 

232 

232540 

5918.3 

-0.19 

23S 

232487 

5917.0 

0.22 

\ 

238 

232439 

5915.8 

-0.09 

2S1 

232163 

5908.7 

-0.15 

''.••I 

257 

232065 

5906.2 

-0.32 

*  '  ” 

260 

231995 

5904.4 

0.05 

.  .i 

Shear 

211 

127776 

3252.0 

0.38 

ft 

(Par.) 

218 

127680 

3249.6 

0.11 

'  J 

221 

127628 

3248.2 

0.03 

*  -  •  ’  ^ 

223 

127599 

3247. S 

-0.29 

■  ‘  ■  1* 

226 

127561 

3246.5 

-0.31 

;  '.1 

234 

127432 

3243.3 

-0.24 

; 

243 

127324 

3240.5 

-0.20 

2S1 

127151 

3236. S 

0.04 

260 

127003 

3232.3 

0.05 

l 

262 

126957 

3231.2 

0.12 

Shear 

210 

127853 

3254.0 

0.09 

(Perp.) 

214 

127790 

3252.4 

0.04 

/•  ■ 

221 

127668 

3249.3 

-0.05 

225 

127611 

3247.8 

-0.28 

230 

127523 

324S.6 

-0.26 

*- 

233 

127473 

3245.3 

-0.35 

236 

127415 

3242.8 

-0.38 

P-  ,  ~v.'< 

241 

127317 

3240.3 

-0.12 

246 

12724S 

3238. S 

-0.04 

•_*  » 

249 

127187 

3237.0 

-0.08 

►  -  • 

rzz 

This  accuracy  in  the  temperature  dependence 
measurements  agrees  very  well  with  the  accuracy 
of  il.6\  determined  by  Barberl<>  in  his  study  on 
aluminum. 

It  can  also  be  noted  from  Table  (2)  that  at 
zero  applied  stress,  the  values  of  dV/dT  for 
shear-parallel  and  shear-perpendicular  are  equal 
within  the  accuracy  of  measurements.  This  means 
that  at  zero  applied  stress,  dV/dT  of  shear 
waves  in  the  steel  investigated  is  texture 
insensitive. 

CALCULATION  OF  STRESS  DISTRIBUTIONS  IN 
SPECIMEN  -  Figure  9  shows  the  calculated  distri¬ 
butions  of  the  circumferential,  the  radial  and 
the  axial  stress  variations  along  a  radial  line 
in  the  specimen  used  in  this  investigation.  The 
stress  distributions  are  calculated  using  the 


relationships. 
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where  o,  is  the  tangential  stress,  or  is  the 
radial  stress  and  ofl  is  the  axial  stress.  U  is 
the  Young's  modulus,  6  is  the  extent  of  inter¬ 
ference,  v  is  the  Poisson's  ratio,  b  is  the  outer 
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l:ig.  S  Typical  variations  of  ultrasonic  longitudinal,  shear  (parallel)  and 
shear  (perpendicular)  velocities  with  temperature  at  0.64  ca.  from 
the  hole  in  reference  speciaen 
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A533B  STEEL 

;:v- 

SPECIMEN  •  1 

’ ."  % 

■ 

B 

LONG.  VELOCITY 

***** 

Variations  of  temperature  dependence  of  ultrasonic  longitudinal  velocity 
as  a  function  of  distance  froa  the  hole  in  speciaen  II.  The  propagation 
is  in  axial  direction.  The  variations  are  shown  before  and  after  plugging 
the  hole. 


Fig.  7  Variation  of  tenperature  dependence  of  ultrasonic  shear  velocity  as  a  function 
of  distance  fro*  hole  in  speciaen  *1.  The  direction  of  polarization  is 
parallel  to  the  reference  line  and  the  direction  of  propagation  in  the  axial 
direction  (shear-parallel).  The  variations  are  shown  before  and  after 
plugging  the  hole. 
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Table (2)  Variations  in  (dV/dT)  obtained  on  reference  specimen  using 
longitudinal,  shear  (parallel),  and  shear  (perpendicular) 
waves . 


Date 

Condition 
of  spec. 

Type  of 

Nave 

dV/dT 

\  Dev 

4/20/82 

As  reed. 

Longitudinal 

-0.640 

1.84 

4/03/82 

As  reed. 

Longitudinal 

-0.652 

0.0 

7/02/82 

As  reed. 

Ixingidutinol 

-0.613 

0.0 

9/23/82 

As  reed. 

Longitudinal 

-0.629 

2.54 

6/05/82 

Machined 

Shear  (Par.) 

-0.401 

0.0 

6/18/82 

Machined 

Shear  (Per.) 

-0.407 

1.47 

7/03/82 

Machined 

Shear  (Par.) 

-0.474 

0.0 

7/14/82 

Machined 

Shear  (Per.) 

-0.483 

0.86 

12/30/82 

Hole  Made 

Shear  (Par.) 

-0.443 

0.0 

12/22/82 

Hole  Made 

Shear  (Per.) 

-0.445 

0.44 

Std.  Dev.  *  0.0172 
»  1.7% 


Theoretical  variations  of  tangential, 
radial  and  axial  stresses  as  a  function 
of  distance  from  the  hole  in  specimen 
•  1. 


diaaeter  of  the  plug  and  c  is  the  outer  diaaeter 
of  the  disc.  The  distributions  shown  in  Fig.  9 
are  calculated  using  the  dimensions  shown  in 
Fig.  I  and  a  radial  interference  of  0.001  in. 
Froa  equations  (5),  (4)  and  (5)  it  can  be  seen 
that  the  shape  of  the  stress  distribution  is  not 
affected  by  the  actual  values  of  interference, 
though  the  magnitude  of  the  stress  components 
depends  on  the  extent  of  that ^quantity . 

TEMPERATURE  DEPENDENCE  OF  LONGITUDINAL 
ULTRASONIC  VELOCITY  -  Figure  6  shows  the  varia¬ 
tions  of  the  temperature  dependence  of  ultrasonic 
longitudinal  velocity,  as  a  function  of  radial 
distance  froa  the  hole  before  and  after  specimen 
was  plugged.  In  these  measurements  longitudinal 
waves  are  propagated  along  the  axial  direction  of 
the  specimen.  From  this  figure  it  can  be  noted 
that  both  before  and  after  plugging,  dV/dT  re¬ 
mains  constant  along  the  radial  direction  of  the 
specimen.  Also  froa  the  figure,  it  is  seen  that 
(dV/dT)  remains  constant  at  values  between 
-0.64  m/s.k  and  -0.63  a/s.k.  This  value  is  very 
close  to  the  value  of  -0.645  obtained  by  Salnma 
and  Nang11  on  a  stress-free  specimen  of 
AS33B  steel.  This  equality  indicates  that  the 
amount  of  the  axial  component  of  residual  stress 
in  specimen  in  asreceived  condition  should  be 
very  saa 1 1 . 

Furthermore,  one  can  sec  from  Fig.  6  that 
after  stresses  are  introduced,  the  axial  com¬ 
ponent  of  stress  also  remains  constant  along 
the  radial  distance  from  the  hole.  This  qualita¬ 
tively  agrees  with  eqn.  (5)  where  og  does  not 
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depend  on  radial  distance  froa  the  hole.  It  is 
difficult,  however,  to  obtain  quantitative 
values  for  the  corresponding  stress  coaponent 
since  (dV/dT)  aeasureaents  are  affected  by  the 
axial  coaponent  of  the  stress,  which  is  parallel 
to  the  direction  of  wave  propagation,  and  cali¬ 
bration  curves  available  are  only  for  the  case 
of  stress  applied  in  a  direction  perpendicular 
wave  propagation. 

TEMPERATURE  DEPENDENCE  OF  ULTRASONIC 

VELOCITY  SHEAR- PARALLEL  AND  SIILAK- 

PERPEND1CULAR  -  It  can  be  seen  from 
equations  (4)  and  (S)  and  Fig.  9  that  the  radial 
and  circumferential  stresses  vary  inversely 
with  the  square  of  the  radial  distance  froa  the 
hole.  Figure  7  shows  the  relationship  between 
the  temperature  dependence  of  ultrasonic  shear 
velocity  (dV/dT) • i  and  the  radial  distance, 
when  the  polarization  is  parallel  to  the 
radial  line.  From  this  figure,  one  can  see 
that  the  data  points  for  the  distances  ranging 
from  1.27  cm  to  3.74  cm  may  be  represented  by 
an  inverse  square  relationship.  However,  the 
data  points  less  than  1.27  cm  show  a  decrease 
in  dV/dT  as  the  distance  decreases.  This  de¬ 
crease  might  be  due  to  the  plastic  deformation 
in  that  region  resulted  froa  the  introduction 
of  the  plug  into  the  hole.  Similar  behavior 
between  (dV/dT),  and  the  radial  distance,  is 
obtained  when  the  polarization  is  aligned 
perpendicular  to  the  radial  line.  These  re¬ 
sults  are  shown  in  Fig.  8.  In  this  figure 
the  behavior  of  the  temperature  dependence  of 
ultrasonic  shear  velocity  is  also  similar  to 
that  observed  when  the  polarization  is  parallel 
to  the  reference  line. 

Again,  it  is  difficult  to  convert  these 
results  of  (dV/dT)g,|  to  stress  quantities  of 
tangential  and  radial  components  of  stress. 

The  measurements  of  dV/dT  in  the  radial  and 
the  tangential  directions  will  be  affected  by 
the  component  of  stress  parallel  to  propaga¬ 
tion,  in  addition  to  that  of  the  component 
perpendicular  to  propagation.  Since  the  cali¬ 
bration  curves  available  are  only  for  the  case 
of  stress  applied  in  a  direction  perpendicular 
to  the  wave  propagation,  it  will  be  difficult 
to  quantitatively  evaluate  the  corresponding 
values  of  stress  from  the  values  of  (dV/dT)  gj. 
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ABSTRACT 

The  temperature  dependence  of  longitudinal  ultrasonic  velocity 
in  polycrystalline  cubic  structures  has  been  derived  in  terms  of 
second-,  third-  and  fourth-order  elastic  constants.  The  derivation 
is  made  using  the  Voigt  average  of  elastic  constants  along  with  the 
expressions  obtained  by  Hiki  et  al.  for  the  temperature  dependence 
of  the  isothermal  second-order  elastic  constants  of  single  crystals. 
Comparisons  between  derived  temperature  dependence  and  measured 
quantities  in  aluminum  and  copper,  yield  values  of  fourth-order 
elastic  constants  which  agree  with  those  previously  obtained  by  other 
investigators  using  single  crystal  data.  Also,  the  comparison 
between  the  calculated  stress  derivative  of  the  temperature 
dependence  and  the  corresponding  experimental  values,  yields 
reasonable  values  for  the  stress  dependences  of  the  third-  and  the 
fourth-order  elastic  constants  of  aluminum  and  copper.  Furthermore, 
based  on  the  fact  that  these  values  are  constants  for  each  material, 
it  is  concluded  that  the  temperature  dependence  of  ultrasonic 
velocity  in  these  materials  varies  linearly  with  stress  as  it  is 
observed  experimentally. 


INTRODUCTION 


A  great  deal  of  experimental  work  has  been  done  on  the 

temperature  dependence  of  second-order  elastic  constants  of  single- 
and  poly-crystals. W  The  data  show  that  the  elastic  constants  are 
linear  functions  of  temperature  up  to  melting  temperatures^.  Also, 
recent  experiments  performed  on  aluminum9>10,  copperll  and  steel12-l4 
elastically  deformed,  show  that  the  ultrasonic  velocity,  in  the 
vicinity  of  room  temperature  changes  linearly  with  temperature,  and 
the  slope  of  the  linear  relationship  is  strongly  influenced  by  the 
presence  of  stress.  In  these  materials,  the  relative  changes  of  the 
temperature  dependence  of  longitudinal  velocity  are  found  to  increase 
by  231  at  a  stress  of  96  MPa  in  aluminum,  by  61  at  a  stress  of  240 

MPa  in  copper,  and  by  a  15*  at  a  stress  of  240  MPa  in  type  A533 

steel.  The  results  on  these  materials  also  indicate  that  the 
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*  relative  change  In  the  temperature  dependence  is  a  linear  function  of 
applied  stress,  and  the  slope  of  this  linear  relationship  appears  to 
remain  unchanged  for  specimens  of  the  same  material  tested.  No 
theoretical  explanation,  however,  is  yet  available  for  this  linear 
behavior  of  the  temperature  dependence  of  ultrasonic  velocity  with 
stress. 


THEORY 


Basically  the  temperature  dependences  of  elastic  constants  are 
due  to  the  anharmonic  nature  of  the  crystal  lattice  and  are  directly 
related  to  the  higher-order  elastic  constants^,  16.  Hiki,  Thomas  and 
Granatol7  have  derived  expressions  for  the  temperature  dependences  of 
the  isothermal  second-order  elastic  constants  using  a  quasiharmonic 
an i strop ic-continuum  model.  The  theory  takes  account  of  polarization 
and  orientation  dependence  of  phonon  frequencies  and  their  strain 
derivatives,  but  assumes  the  derivatives  are  wavelength  independent. 
The  explicit  expanded  forms  for  cubic  crystals  are  given  as 
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where  the  Gruneisen  parameters  are  given  as. 


(3) 


Y-1  ■  -(l/2w){jjN2+cJ2(H|*H|)+2wuJ+CinN2u2 

+  C144(N2U3+N3U2)2+C112fN2U^3^2NlUl(N2U2^3S^ 

+  2C123N2N3U2U3+C  166t(NlU3+N3Ul)2+(NlU2+N2Ul)2]  J  *  (4) 

Y-3  =  -(1/w){cJ4N2N3+wU2U3^C144N1U1(N2U3+N3U2) 

♦  C166[(N2  +  N3)U2U3  +  N2N3(U2  *  U3)] 
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with 

w  *  c|1(n|u|+n|u|^n|u|)  «•  c|4[(N2u3^N3u2)2  ♦  (n3u1+n1u3)2 
*  (N1U2+N2U1)2]  ♦  2C^(N2N3U2U3*N3N1U3U1^1N2U1U2). 
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In  these  expressions  N  and  I)  are,  respectively,  the  propagation  and 
the  polarization  vectors  of  the  ith  normal  mode,  the  superscripts  S 
and  T  are,  respectively,  for  the  adiabatic  and  the  isothermal  elastic 
constants,  k  is  the  Boltzmann  constant,  and  p  is  the  density  of  the 
undeformed  material.  In  evaluating  the  lattice  sums  of  the  above 
equations,  Hiki  et  a  1 . 1 7  used  a  367-point  grid  over  an  octant  of  the 
Debye  sphere.  The  values  of  the  sums  obtained  by  this  method  for  Cu, 
Ag,  and  Au  have  been  shown,  however,  to  be  only  a  fewer  percent 
larger  than  those  calculated  using  a  pure  mode  average^.  Employing 


the  latter  method  and  neglecting  the  differences  between  C 
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which  have  been  shown  to  be  small,  the  temperature  dependences  of  the 


elastic  constants,  Cllf  c*  “-^(Cn  -  C12)  and  can  be 
approximately  expressed  as 
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Equation  (7)  was  obtained  by  considering  the  longitudinal  waves 
propagating  along  the  1 1 00]  direction,  while  in  evaluating  the 
temperature  dependences  3C/0T  and  **C../3T,  we  considered  the  inodes 
of  shear  waves  propagating  along  theulO]  direction  and  polarized 
parallel  to  the  [110]  and  [001]  directions,  respectively. 

The  expressions  for  the  temperature  dependences  of  the  second- 
order  elastic  constants  show  a  linear  relationship  up  to  the  fourth- 
order  elastic  constants.  Assuming  the  Cauchy  relations??  hold  for 
the  fourth-order  elastic  constants  as, 


C1 112  =  C1 155  ’  C1122  *  C1266  “  C4444* 

CH23  *  C1144  “  C1255  *  C1456  *  C4455 


then  only  four  of  the  fourth-order  elastic  constants  namely,  Cim, 


C-| 1 12 *  Ci  1 22  and  Cii23  need  to  be  considered.  Basically  the  Cauchy 
assumption  means  that  the  short-range  repulsive  forces,  which 
contribute  the  most  to  the  higher-order  elastic  constants,  may  be 
reasonably  represented  by  central  forces.  Also,  assuming  that  only 
nearest-neighbor  interactions  contribute  to  the  fourth-order  elastic 

constants20,  we  obtain 


C1112  *  C1122  *  TC1111  and  C1123  * 


and  for  third-order  elastic  constants, 


Clll  ■  2C112  "  2C166 
C123  *  C4S6  *  C144  '  0 


This  assumption  follows  from  the  predominant  contribution  of  the 
short-range  forces  to  these  higher-order  elastic  constants. 

Accordingly  eqns.  7,  8,  and  9  can  respectively  be  written  as 


3C11 

<TT>  * 


-kpo„3Cll  4  Clll,2  5C111  4  Cllll, 

9  l  l  r  /  p  J 

li  Ln 


^ 3T  J  *  ^ 


111  4  cllll] 


3C.. 

( — 51) 

1  3T  J 


'111  4  Cllll^ 


Garbar  and  Granato?1  used  these  relationships  to  determine  values  of 
the  fourth-order  elastic  constants  of  noble  metals  Cu,  Ag  and  Au. 
The  values  of  the  fourth-order  elastic  constants  obtained  were 
approximately  an  order  of  magnitude  greater  than  and  opposite  in  sign 
from  Ciii,  which  agrees  with  the  predictions  of  nearest 
interact ions22  for  a  potential  of  the  form  1/rn  where  Cim/Cm  ■ 

■'“(n  +  6)  and  a  value  of  13  was  suggested  for  n. 


TEMPERATURE  DEPENDENCE  OF  POLYCRYSTALS 


Voigt23  has  shown  that  the  Young's  modulus  and  the  rigidity 
modulus  of  a  polycrystal  may  be  calculated  from  the  single-crystal 
elastic  moduli  if  one  assumes  that  each  grain  is  deformed  to  the  same 
strain.  His  results  yield, 


<C11-C12  *3C44Hc„*3C12) 

3Cj2  ♦  C44) 


(15) 


C11  *  C12  +  3C44 
5 


(16) 


The  velocity  of  longitudinal  waves  in  a  polycrystal  is  related  to  the 
Young's  modulus  and  the  modulus  of  rigidity  as 


1/2 


(17) 


where  p  is  the  density  of  the  solid. 

Substituting  eqns.  15  and  16  into  eqn.  (17),  yields 


V£  = 


/5p 


(5Cn  -  4C'  ♦  4C^] 


1/2 


44J 


(18) 


where  C*  is  Cjj  -  C12/2. 

Considering  the  density  p  remains  unchanged  over  a  short  period 
of  temperature  and  the  change  in  the  ultrasonic  velocity  with 
temperature  is  due  only  to  changes  in  the  elastic  constants,  the 
temperature  dependence  of  ultrasonic  velocity  will  be  given  by, 


2/§p" 


-1/2  ^11 
(5Cn-4C’*4C44)  1,1 [S-yf- 


(19) 


Substituting  for  the  temperature  dependences  3C../3T,  3C'/3T  and 
dCj^/dT  from  eqns.  (12,  13  and  1h),  eqn.  (19)  takes  tne  form, 


^(SCn-4C-.4C44,-^,-^-^ 


'‘'Ill 

4C,, 


‘'"'111 

16C' 


Equation  (20)  represents  the  temperature  dependence  of  longitudinal 
velocity  in  a  polycrystalline  solid  in  terms  of  the  second-,  third- 
and  fourth-order  elastic  constants  of  the  single  crystal  structure  of 
that  solid.  A  similar  relationship,  however  more  complicated,  may 
also  be  obtained  using  the  Reuss24  averaging  procedure  in  calculating 
the  Young's  and  Rigidity  moduli,  which  assumes  each  grain  is  subject 
to  constant  stress. 


CHANGE  OF  TEMPERATURE  DEPENDENCE  OF  ULTRASONIC  VELOCITY  WITH  STRESS 
1.  Aluminum 

Measurements  of  the  ultrasonic  longitudinal  velocity  in  the 
temperature  range  between  1 80  and  260  K  on  the  aluminum  alloys  2024-0 
and  6063-T4,  have  shown  that  the  velocity -decreases  1  (nearly  -  with 
temperature,  and  the  value  of  the  temperature  dependence  dV  /dT  is  in 
the  range  between  -92.3  and  -111.1  cm/s.k.  Similar  results^were  also 
obtained  on  the  aluminum  alloys  2024-T351,  3003-T251,  and  6O6I7T4, 
where  the  values  of  dV  /dT  are  respectively  found  to  be  -118.7,  - 
125.1  and  -130.4  cm/srk.  The  average  value  of  the  temperature 
dependence  of  longitudinal  velocity  in  polycrystalline  aluminum,  -115 
m/s.k.,  and  the  values  of  the  second-  and  third-order  elastic 
constants  obtained  by  Thomas25,  C  -  1.0675»  C'  -  0.237,  Cw  - 
0.2834  and  Cm  “  “10. 76  in  units  ot  1012  dyne/cm  ,  are  substituted 
into  eqn.  (20).  The  value  of  the  fourth-order  elastic  constant  Cmi 
is  found  to  be  +86.5  x  1012  dyne/cm  ,  which  is  opposite  in  sign,  and 
about  a  factor  of  8  larger  than  that  of  measured  Cjjj.  This  agrees 
with  the  predictions  of  nearest-neighbor  interactions  for  a  potential 
of  the  form  1/rn  with  n  ■  10.  It  also  agrees  with  the  findings  of 
Garbor  and  Granato  in  their  determination  of  the  fourth-order  elastic 
sonstants  of  the  noble  metals  Cu,  Ag  and  Au. 

Equation  (20)  indicates  that  the  temperature  dependence  of 
ultrasonic  velocity  is  a  constant  quantity  which  may  be  computed 
knowing  the  values  of  the  second-,  third-  and  fourth-order  elastic 
constants.  Again  considering  the  density  remains  unchanged  as  a 
function  of  stress,  the  derivative  of  the  right-hand  side  of  eqn. 
(20)  with  respect  to  stress  will  be  a  function  of  the  second-, 
third-  and  fourth-order  elastic  constants,  and  their  derivatives  with 
respect  to  stress,  3Zvo/jT3o  will  then  be  a  constant  as  long  as  the 
derivatives  of  the  elastic  constants  with  respect  to  stress  are 
constants.  Recent  experiments  on  aluminum,  copper  and  steel,  have 
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shown  that  the  temperature  dependence  of  ultrasonic  velocity  In  these 
materials  is  a  linear  function  of  applied  stress.  The  slope  of  this 
linear  relationship,  however,  is  to  be  different  in  different 
materials,  and  differs  according  to  the  relative  direction  of  stress 
with  respect  to  the  directions  of  propagation  and  polarization. 

.  In  order  to  obtain  quantitative  determination  of  the  derivative 
®  vfc/dTdo  a  computer  program  was  developed  to  compute  the  value  of 
the  quantity  ®^v  /3o3t  for  given  values  of  second-,  third-and 
fourth-order  elastic  constants  and  their  derivatives  with  respect  to 
stress.  The  program  is  based  on  Newton's  backward  difference  formula 
for  polynomial  approximat  ion26f  which  does  not  give  the  exact 
differential  but  rather  an  approximate  differentiation  of  the 
function.  The  accuracy  of  the  program  was  tested  for  functions  of 
known  derivatives,  and  the  results  of  the  derivatives  at  given  values 
of  the  function  were  found  to  be  equal  to  those  of  the  exact  solution 
within  0.1S. 


The  second-  and  the  third-order  elastic  constants  of  aluminum 
single  crystals  obtained  by  Thomas  along  with  the  value  of  the 
fourth-order  elastic  constants,  +86.5  x  1012  dyne/cnr  obtained  from 
t^e  temperature  dependence  calculations,  are  used  to  evaluate 
3  v£/3o3j.  Also  the  values  of  the  stress  dependences  of  the  second- 
order  elastic  constants  Sc^/Bo,  3c*/3o  and  for  stress 

applied  perpendicular  to  propagation,  obtained  by  momas25  are  used 
in  the  stress  computer  program  to  calculate  3zV£,/3a3T.  These  values 
for  stress  applied  in  tension  are  +1.67,  -0.61  and  -2.46 

respectively.  The  value  —0.6 1  is  the  average  of  -0.32  and  -0.90 
obtained  with  stress  applied  in  two  different  directions.  For  stress 
dependences  of  the  third-  and  fourth-order  elastic  constants  3C^i/3o 
and  ®Cjjii/3°,  no  values  are  available  in  the  literature,  and 
accordingly  some  estimates  were  made.  Using  3C|^/3o  -  10  and 

^^111l/^a  *  J00  is  .found  to  yield  a  value  of  0.228  x  10”7 
cm/s.k/dyne/cm  for  3zv8,/3o3t  which  is  .an  excellent  agreement  with 
the  value  0.242  x  10“7  cm/s.k/dyne/cm^  obtained  experimental  1 yll. 
The  values  of  ^Cm/3o  and  ^Cmi/So  used  to  yield  the  agreement 
between  calculations  and  experiment  are  of  the  same  order  of 
magnitude  to.  be  expected  for  these  quantities.  The  values  of  the 
derivative  3zv8,/3o3j  »  0.12  x  10"7  cm/s.k/dyne/cnr  when  both  3Cm/3o 
and  ^Cuji/So  are  equal  to  zero.  Also  changing  the  values  of 
®Cin/3o  and  3cmi/do  by  *5Wi  changes  the  value  of  the  derivative  by 
*25$.  Table  I  contains  the  values  of  §4^/3o3t  at  different  values 
of  ®C|ii/3o  and  3Z\A/3o3t  at  different  values  of  3Cm/3o  and 
3Cn,i/3a. 


Table  I. 


Values  of  the  derivatives  d  v&/dodT  in  aluminum  at 
various  values  of  stress  dependences  C^/do  and 
dCuu/do.  Calculations  are  made  using  Cn  -  1.068, 
C'  -  0.2317,  «  0.2834,  Ciii  x  -10.76  and  Cim  “ 

86.5  in  units  of  1012  dynes/cm  ,  dCjj  /Bo  ■  I.67, 
ac/ao  -  0.61,  ac^/aa  -  -2.46. 


3C111 

aCllll 

dv 

3  ,  £ 

^  cm/s.k 

do 

do 

do  C  dT 

1  2 
dynes/cm 

0 

0 

0.116  x 

io-7 

s 

50 

0.172  x 

ID'7 

10 

100 

0.228  x 

io-7 

15 

150 

0.284  x 

io-7 

Experiment 

0.242  x 

io'7 

2.  Copper 

The  second-  and  the  •  third-order  elastic  constants  of  single 
crystal  copper  have  been  measured  by  Hik i  and  Granato20.  The  values 
obtained  are  C^i  «  1.661,  £'  ■  0.231,  C.j,  *  0.75&  and  Cm  -  -12.71 
in  units  of  1 0 »2  dynes/c<n  .  Also  the*H temperature  dependence  of 
ultrasonic  longitudinal  velocity  3v£/3T  in  CDA  110  polycrystalline 
copper  has  been  measured  by  Salama  and  Lingll.  The  results  obtained 
on  annealed  specimens  averaged  a  value  of  -0.4865  m/s.k  for  (3v£/3t), 
while  those  on  as  received  specimens  have  an  average  of  -0.521  m/s.k. 
Substituting  the  average  of  these  two  quantities  -0.504  m/s.k  for 
3v£/3t  in  eqn.  (20),  and  using  the  values  of  the  second-  and  third- 
order  elastic  constants  obtained  by  Hlki  and  Granato20,  we  obtain  a 
value  of  C|]]|  ■  93.55  x  10l2dyne/cm.  This  value  is  opposite  in 
sign  and  about  a  factor  of  7*5  larger  than  that  of  the  third-order 
elastic  constant  Cju.  It  also  lies  within  the  range  87.7  ~  107  x 
10l2dyne/cmz  obtained  by  Garbar  and  Granato21  from  their  analysis  of 
the  temperature  dependence  of  the  second-order  elastic  constants  of 
single  crystals  of  copper. 

The  values  of  the  stress  dependence  of  the  second-order  elastic 
constants  Cn,  C'  and  Crr  in  copper  were  calculated  using  the 
measurements  of  Hiki  and  Granato20f  and  found  to  be  dCjj/ds  *  -0.66, 
dC'/dcJ  -  -1.17  and  dC../do  -  1.41  for  tensile  stress.  Each  of  the 
values  of  dC'/do  and  dCj./do  are  the  average  of  the  values  obtained 
on  three  different  combinations  of  stress,  propagation  8nd 
polarization  directions.  Each  of  these  two  quantities  depend  strongly 
on  the  relative  orientation  of  the  direction  in  which  the  uniaxial 
stress  is  applied  with  respect  to  the  polarization  direction.  Using 


these  values  along  with  those  of  the  second-,  third  and  fourth-order 
elastic  constants  of  copper  in  the  computer  program  to  obtain  the 
derivative  3  a  good  agreement  between  theory  and  experiment 

for  this  quantity  is  obtained  with  C m/*0  -  8  and  aCmi/*°  ■  80. 
The  calculated  value  of  3z\£^o3T  is  found  to  be  0.129  x  10"® 
cm/s.k/dyne/cm  ,  while  that  obtained  from  measurements  of  the 
temperature  dependence  of  ultrasonic  velocity  as  a  function  of 
applied  stress  is  equal  to  0.125  x  10“8  cm/s.k/dyne/cm2.  These 
values  are  a  factor  of  20  smaller  than  those  obtained  in  aluminum. 
Nevertheless,  the  values  of  the  stress  dependences  of  the  third-  and 
fourth-order  elastic  constants  3  Cm/00  -  8  and  3Cim/3o  «  80 

obtained  in  copper  are  of  the  same  magnitude  as  of  10  and  100 
computed  for  °Cm/^a  and  3Ciiil/®°  respectively  in  aluminum. 
Contrary  to  aluminum,  however,  the  calculated  value  of  the  slope  of 
the  relationship  between  the  temperature  dependence  of  the 
longitudinal  velocity  3%1/9t  and  applied  tensile  stress  o,  is  found 
to  be  sensitive  to  the  values  used  for  3cm/3 a  and  3Cmi^°-  Table 
II  lists  the  calculated  values  of  32v&/3o3t  in  copper  using  different 
values  of  3c^i/3o  and  3Cmi^°  along  with  the  same  values  of  the 
second-,  third-  and  fourth-order  elastic  constants  and  the  stress 
derivatives  of  the  second-order  elastic  constants. 


Table  II 


Values  of  the  derivative  3zvJ,^t3o  in  copper  at 
various  values  of  the  stress  dependences  3Cm/^°  and 
**c111l/d°«  Calculations  are  made  using  Cf i  *  1.661, 
C'  -  0-231,  C.j.  -  0.756,  Cm  s  -12.71  and  Clni  “ 
93-55  in  units  of  10^2  dyne/cm.  3Cji/3 a  ■  -0.66, 
3c'/3o  .  -1.17  and  3C../3 a  -  i.4l. 
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